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ABSTRACT 

The characteristics of light variation of RSGs in SMC are analyzed based on 
the nearly 8-10 year long data collected by the ASAS and MACHO projects. 
The identified 126 RSGs are classified into five categories accordingly: 20 with 
poor photometry, 55 with no reliable period, 6 with semi-regular variation, 15 
with Long Secondary Period (LSP) and distinguishable short period and 30 with 
only LSP. For the semi-regular variables and the LSP variables with distinguish- 
able short period, the K s band period-luminosity (P-L) relation is analyzed and 
compared with that of the Galaxy, LMC and M33. It is found that the RSGs 
in these galaxies obey similar P-L relation except the Galaxy. In addition, the 
P-L relations in the infrared bands, namely the 2MASS JHK S , Spitzer/IRAC 
and Spitzer/MIPS 24 fj,m bands, are derived with high reliability. The best P-L 
relation occurs in the Spitzer/IRAC [3.6] and [4.5] bands. Based on the compar- 
ison with the theoretical calculation of the P-L relation, the mode of pulsation 
of RSGs in SMC is suggested to be the first overtone radial mode. 

Subject headings: stars: late-type — stars: oscillations — stars: variables: other — 
supergiants 
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Introduction 



Red Supergiant stars (RSGs) are massive evolved, He-burning, extreme Population 
I stars with mass range roughly between 10-30 M Q . They have relative low effective 
temperatures of about 3450 ~ 4100 K (spectral type of late-K ~ M) , enormous radii 



up to 1500 Rp and high luminosit ies of about 2000 ~ 30000 L (ILevesque et al. 



Massey et al. 



2005 



2008 



Levesque 



20101 ) . The very large radii make them one class of the largest 



stars in the universe in terms of physical size, although not the most luminous or most 
massive stars. Due to the large radii, their gravity is not strong enough to hold the matter 
in the outer layer, which gives rise to large mass loss rate (MLR) in the order of magnitude 
of about 10~ 6 to 10~ 4 M yr _1 . This strong wind cr eates a thick dusty envelope that makes 



it very hard to determine their physical parameters ( 



Groenewegen et al. 



2009 



Yoon fc Cantiello 



2010 



Josselin et al. 



Sargent et al. 



2000 



Massey et al. 



2005 



201 lh . Further more, RSGs 



have long been known for their hundreds-of-days semi-regular light variatio n, as a SRc class 
in General Catalogue of Variable Stars (GCVS; iKholopov et al.lll985-1988l ). which may be 



1994 



( Stothers 


1969; 


Wood. Bessell & Fox 


1983 


: Low et al. 


1984 


Heger et al. 


1997; 


Bono & Panagia 


2000 


; < 


jluo & Li 


2002 


)• I 



Schaller 



1990 



Li fc Gong 



20021 ) . Besides, they also show an 



inexplicable long-term variation with a dura tion up to 4000 days or longer, known as the 



Long Secondary Period (LS 



Messina 



2007 



Percy Sato 



) ph enomenon (IStothers 



1972 



Kiss. Szabo fc Bedding 



2006 



20091 ). Up to now, the mechanism of LSP is still unclear. The 



models of binary, pulsation, convection cell and surface hot spot are st udied but none of 



them fully agrees with all the observat i ons and theoretical 



Wilson et al. 



1992 



Tuthill et al. 



1997 



Groenewegen 



2004; 



e xpectations ( 



Buscher et a 



Wood. Olivier fc Kawaler 



1990 



2004 



Kiss . Szabo fc Bedding 



2006 



Messina 



2007 



Haubois et al. 



2009; 



Nie et al. 



2010 



Stothers 



20101 ). In addition to the semi-regular short-term and LSP variations, RSGs yet presents an 
irregular photometric variation that is a feature of the Lc class in GCVS. It can be explained 
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by the large convective ce 



2006 



simulations (ISchwarzsch ild 



Chiavassa et al 



s in t 



2010 



975 



le surface of the star that are proved by some recent 



Ant ia. Chitre &: Narasimha 



Arnett fc Meakin 



2011) 



1984! : 



Kiss. Szabo fc Bedding 



A period-luminosity (P-L) relation is found in the light variation of RSGs, which 



is important for the potential use o 
of variable s was fir s t discovered by 



applied by 



Hubblel ( 11925 



RSG a s a ve ry luminous candle. The P-L relation 



Leayitt 



(]1908[ ) for the Cepheid variables and then 



19261 ) to demonstrate unambiguously that NGC 6822 and 



M33 were extragalactic systems. Since then, Cepheids have become a major tool to 
determine the extragalactic distance. RSGs, as one class of the most lum inous stars, 
can extend the scale of distance even further than Cepheids. iGlass! (119791 ) discussed the 



probability of usin g luminous RSGs as an extragalactic distance indicator originally. Then, 



Feast et al. 



( 1l980l ) used infrared photometry to derive the bolometric magnitudes and 



demonstrated the existen ce of P-L relation for RSGs in the Large Magellanic Cloud (LMC). 



Catchpole &: Feast 



(119811) confirmed the relat i onship in the Small Magellanic Cloud (SMC) 



in the same way. Later, IWood. Bessell fc Fox! ( 119831 ) made use of infrared JHK photometry 



and low-dispersion red spectrum to investigate the Long Period Variables (LPVs) in the 
LMC and SMC. They found that the LMC RSGs had a similar P-L sequence to the 
Asymptotic Giant Branc h stars (AGBs) but were approx i mately one m a gnitu de br ighter in 



the K band. Afterwards, 



( 11993 



Kinman. Mould fc Wood! (119871 ). 



Mould et al. 



(11990( 1 and 



Dambis 



utilized the infrared P-L re 



ationof RSGs to estimate the distance to the extragalaxy 



M33. iPierce. Jurcevic fc Crabtred (120001 ) re-calibrated the P-L relations of RSGs in the 
Perseus OBI (Per OBI) association, LMC and M33 in various bands and suggested a 
universal relation in these h eterogeneous environments to measure the distance to M101 



resulting in about 6.4 Mpc. 



Kiss. Szabo fc Bedding! ( 120061 ) obtained a P-L relation similar 



to that of AGBs in LMC based on the analysis of 48 Galactic RSGs with almost 60-year 
data. After all, we (Yang & Jiang 2011, Hereafter Paper I) made use of visual time-series 
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photometric data collected by the ASAS and MACHO projects that span nearly 8-10 years, 
and established the P-L relations of RSGs in LMC in the infrared bands including the 
2MASS JHKs bands and Spitzer /IRAC and MIPS 24/xm bands. Our study proved that the 
semi-regular RSGs variables and the LSP variables with distinguishable short period show 
a tight P-L relation in the infrared bands but relatively sparse in the V band which may be 
caused by the inhomogeneous interstellar extinction. The results are compared with other 
P-L relationships for RSGs and the P-L sequences of red giants in LMC. The discrepancy 
between various results may be due to sample selection, the error in the periods or the 
indicator of the luminosity. For further improvement, we extend our work from LMC to 
SMC and try to find a consistent result within our Galaxy and MCs based on similar data, 
similar criteria for selecting RSGs and analyzing methods. 



2. Sample Selection and Data Analysis 



A sample of RSGs pure and as large as possible in SMC is the basis to derive a reliable 
P-L relation. As a first step, we assembled the candidates of RSGs from previous works 
as our targets to consist a large sample. One should be cautious for the selection because 
RSGs are easily confused with AGBs since both are r e d and luminous while RSGs are a 



bit more luminous and less red. 



Wood. Bessell fc Fox 



(11983 1 suggested that M bol = -7.1 



can be a criterion to discriminate AGBs and RSGs. But sometimes the most luminous 



AGBs known as super-AGBs can be so bri ght as Mbol up to —8 mag (ISiess 



Eldridge et al. 



2007 



2006 



2007 



Poelarends et al. 



20081 ) . On the other side, some low- luminosity RSGs 



surrounded by a thick dusty envelope can be darkened to the average luminosity level of 
AGB stars. In addition, previous works selected RSGs in different way s. Early in 1981, 



i nfrare d photometry of 22 RSGs variables in the SMC was performed by 



Catchpole fc Feast 



( Il98ll ). Excluding two anomalous objects, their fitting to 20 targets derived a good P-L 



relation as iFeast et al.l (119801 ) did for LMC. iWood. Bessell fc Foxl (119831 ) produced a new 
catalog of 42 LPVs in SMC from a detailed analysis of their near-infrared photometry and 
low-dispersion red spectra. But according to their criterion for AGBs an d RSGs, there were 



only 5 targets identified as RSGs, which is discrepant with the result of 



Catchpole fc Feast 



( 1l98ll ). A simple and easy solution to this discrepancy is to shift the borderline between 
AGBs and RSGs from Mt, i = —7.1 to Mh n i = —6.1 , whic h would enlarge the sample size 



to 22 objects and coincide with 



sample in SMC was made by 



Catchpole fc Feast 

J 



jl98lh . A big extension of the RSGs 



Massey fc Olsenl (120031 ) who selected 105 s ources through 



multi-object spectroscopy of a sample of red stars pr eviously ident ified by 



separating bona fide RSGs from foreground dwarfs (IMassey 



Massey 



(I2002h in 



19981 ). They simultaneously 



made use of high-accuracy (< lkms l ) r adial velocities f or all the c andidates to confirm 



them as RSGs. Also based on this work, 



Levesque et al. 



(12006f l and 



Massey et al. 



(120071 ) 



identified additional 9 sources. The combination of the three s amples consists the mo st 



van Loon et al 



(120081 ) 



comprehensive and reliable sample of RSGs in SMC. Recently, 
identified 8 RSGs in SMC acoording to their K-band luminosity and the VLT/ISAAC 3-4 
fim spectra th at showed the first oy ertone SiO band at 4 /im in emission (see Fig. 6 in their 



paper). They (Ivan Loon et al.l 120101 ) also presented Spitzer far- infrared spectra between 52 



93 /im by MIPS for 2 RSGs to investigate their infrared and dust properties. 



Ou r preliminary sample i s a combination of all th e 



1981) 



Wood. Bessell fc Fox ( 



Massey et al 



(120071 ) 



1983) 



van Loon et al. 



Massev fc Olsen (2003) 



sam pl es from 



Catchpo 



Levesque et al 



van Loon et al. 



e & 



(120061 ). 



cast 



(120101 ). which brings about 



(120081 ) and 

166 objects. Considering the redundancy between some samples, the same procedure as in 
Paper I is taken, i.e. to adopt the sources for its first appearance in paper published and 



subtract previous sources from fo 
of 140 stars with 20 targets from 



lowing papers. Co n seque ntly, the final sample consists 



(119831 ). 101 from 



Catchpole &: Feast 



(1198 if). 6 from 



Massey fc Olsenl (120031 ) , 8 from lLevesque et al. 



Wood. Bessell fc Fox 



(120061 ). 1 from 



Massey et al. 



( 120071 ). 3 from Ivan Loon et al.l (120081 ) and 1 from Ivan Loon et al.l (120101 ) . However, the 
sample will be further purified by the brightness and color indexes later. 



The photometric data i n visual bands ar e taken from the databases of the All Sky 



Automated Survey (ASAS) dPoimans 



(MACHO) projects (lAlcock et al. 



vl 



20021 ) and the MAssive Compact Halo Objects 



19971 ) . The one-epoch near- and mid-infrared photometric 



data ar e retrieved from the T wo Micron All Sky Survey (2MASS) Point Source Catalog 



(PSC) (ISkrutskie et al. 



20061 ) , and the Spitzer/ 'SAGE Legacy Program (IMeixner et al. 



20061 ) database. 



2.1. Color-Magnitude Diagrams and Two-Color Diagrams 

2.1.1. Near- Infrared Color- Magnitude Diagram J — K$ vs. K$ 

As we mentioned earlier, AGBs and RSGs are easily confused if judged only by their 
brightness. For more accurate identification of RSGs, effective temperature has to be 
taken into account. With both brightness and effective temperature, the color-magnitude 
diagrams (CMDs) and two-color diagrams (TCD) are powerful tools to identify RSGs. 
However, RSGs sometimes have thick dusty envelop due to large MLR whose extinction 
affects the position of RSGs in these diagrams. In order to avoid this extinction effect 
as much as possible, the near- and mid-infrared bands of 2MASS, Spitzer /IRAC and 
MIPS are chosen because the extinction in infrared is much smaller than in visual. The 
infrared magnitudes of these candidates are taken from the SAGE-SMC IRAC Combined 
Epoch Catalog vl.5 that already combined the Epoch 1 and Epoch 2 (the SAGE project 
is composed of two-epoch observations separated by three months) IRAC images and 
cross-associated with the 2MASS PSC. For the MIPS 24 /xm band, there is not yet any 
combined catalog released, the SAGE-SMC MIPS 24 /iin Epoch 1 Catalog is then adopted 
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and the missing targets are extracted from the Epoch 2 and Epoch (the IRAC and 
xom t he S 3 MC pathfinder survey of the inner 3 square degree of the 



MIPS observations 


:rom t 


SMC, 


Bolatto et al. 


2007) 



20071 ). These two catalogs are cross-associated and with high-quality 



photometry as the selected subsets of the I RAC and MIPS Arch ive Catalog. More details 



Meixner et al 



(120061 ) and the SAGE Data 



about the SAGE catalogs can be found in 
Description Document^. 

The procedure of identifying the optical RSGs in the IR catalogs is the same as done in 
Paper I. A search radius of 1" is used for the SAGE catalogs that agrees with the nominal 
pointing accuracy of Spitzer. The closest and brightest counterpart is identified from all the 
objects within the search circle. The sources that have null results in the 2MASS or IRAC 
or MIPS bands are retained instead of being removed. Fig. [1] shows the spatial distribution 
of all the 140 sample stars superposed on the Spitzer /SAGE 8 fim mosaic image. It can 
be seen that the sources are mostly located in the 'bar', probably because that the SMC 
'bar' contains more star-forming regions and young stars than the 'wing'. Table [1] lists 
the coordinates and infrared brightness of all the sample targets, where the ' — ' symbol 
means the data is lacking in the corresponding band. In total, 38 measurements are lacking 
for 13 objects, i.e. No. 34 and 132 without infrared data, No. 61 with only the [24] band 
data, No. 15 and 73 with no JHKs band data, No. 40 without the [4.5] and [8.0] band data, 
No.2, 27 and 55 without the [4.5] band data, No. 7, 19, 86 and 108 without the [24] band 
data. These targets will be absent in further analysis of CMD and TCD when related. The 
RSGs candidates are compared with the 3654 massive stars (M ^ 8 M Q ) in the SMC from 
Bonanos et al.l (120101 ) whose infrared properties have been carefully investigated. They 
are collected from literatures and have been identified by the same criteria in the SAGE 
infrared data as ours. 



^ttp:/ /irsa.ipac. caltech.edu/data/SPITZER/SAGE/doc/ 
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The brightness is a key criterion in identifying RSGs. In Paper I, the bolometric 
magnitude is converted fro m the Ks band magni tude by a simple shift of 3 magnitudes as 



m bo i = m K + 3 suggested by 



Josselin et al. 



( 120001 ) (the difference of the brightness between 



the K and Ks bands is very small and ignored in this work). This conversion was used 
for Galactic RSGs but reliable only under specific condition, i.e. the circumstellar and 
interstellar extinctions are already corrected. Considering that the metallicity of SMC 
is much lower than the Galaxy and no correction of the interstellar extinction that is 
comparable to the observational uncertainty (the extinction in the J and Ks band is about 



19911 ) is adopted), a 



0.08 and 0.03 mag respectively if E(B — V) = 0.1 ( ISchwering &: Israel 
different conversion is needed from the Ks band magnitude to the bolometric magnitude, 
especially when there are some heavily dust-enshrounded RSGs. Thus, a new linear relation 
betw een infrared bands and bolometri c magnitude is derived from 113 ta rgets in total, 7 2 



from [Massey &: Olsenl ( 120031 ). 37 from iLevesque et al. 



(]2006[ ) and 4 from ILevesque et al 



( 120071 ) respectively. The bolometric magnitudes of these sources are based on the spectral 
fitting since there may be some large (an d in some cases, unce rtain) reddening if bolometric 



magnitudes are derived from (V — K) (ILevesque et al. 



20061 ). Fig. |2] shows the result 



of fitting at a form of m^ = a x mt,oi + b which indicates a moderate discrepancy from 
nibd = % + 3 shown as a dashed line in the Ks band. Fitting parameters are listed in 
Table [2] which will be used in determining the limits of brightness in the CMDs. 

Hereafter, the ID number in Table [1] is used to code the targets in CMDs and TCDs 
for convenience. Fig. [3] is the near-infrared CMD, the Ks versus J — Ks of all targets. In 
this diagram, most targets are located in a region of 7.5 < Ks < 10 and 0.5 < J — Ks < 1.6, 
which have similar range in J — Ks, but fainter in the Ks band compared to LMC. The 
fainter Ks brightness can be partly explained by the larger distance of SMC, and also 
possibly the fewer number of very massive RSGs in SMC. 
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The boundaries of luminosity and color index are set up in the same way as in 
Paper I with a modification in the boundaries of the Kg magnitude. The new boundaries 
of the Kg magnitudes are 7.6 and 11.0 shown as the dashed horizontal lines in FigJHl 
which is determined based on the mass-lu minosity relation for mas sive stars (9-30 M ), 



L/L = (M/M ) 7 with 7 very close to 4.0 flStothers fc Leung 



197ll ). From the luminosity, 



the bolometric magnitude is calculated by M bo i = 4.74 — 2.51og(L/L ) and the Kg 
magnitude is derived from our newly derived relati onship between niK, and mboi as listed in 



Table [2] with a distance modulus of 18.91 for SMC (IHilditch et al 



20051 ) . The corresponding 



mass of the most targets ranges from 13 M to 30 M , a bit higher than used (10-25 M ) 
in Paper I, which is consistent with the expected mass range of RSGs although lacking 
the low-mass (9-13 M ) RSGs. For the color index J — Kg, the limits of .5 and 1.6 mag 



correspond to the lower observat ional boundary defined by 



boundary of carbon-rich stars by iHughes &: Wood 



Josselin et al. 



(J2000j) and the 



(119901 ) respectively. With the restrictions 



on both the Kg band magnitude and the color index J — Kg, 130 targets (92.9%) are 
qualified. 

Other suggestions on the criteria for RSGs are checked as well. The red limit of 1.2 mag 
in J — Kg is determine d from the evolutionary synthesis models of RSGs in near-infrared 



by 



Origlia et al 



(119991 ) (see Figure 3 and 11 in their paper; although their model fitted 
the metallicity of LMC, there should be no remarkable difference to shift the limit of 1.2 
mag for SMC). The lower limit of Kg = 9.5 mag is suggested to convert from M^i = —7.1 



mag to distinguish the AGBs and RSGs by lWood. Bessell fc Foxl (jl983l ). nevertheless there 
are alternative values of M hni to discriminate RS Gs from AGBs, such as M hni = — 7.5 mag 



by 



van Loon et al. 



Groenewegen et al. 



3). 



(2009). These two 



(119991 ). M bol = -8.0 mag by 
criteria (Kg = 9.5 mag and J — Kg=1.2 mag) are shown in the diagram as dotted lines. The 
modified limit of Kg < 9.5 mag picks out 111 targets (79.3% of the sources), and the limit 
J — Kg < 1.2 mag picks out 126 targets (90% of the sample), both of which are satisfied 
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by most of the objects. There are 103 targets withi n both the Ks m a gnitu de and color 



Massey &: Olsenl (120031 ) catalog that 



index limits, i.e. 73.6% and most of them are in the 
confirmed the reliability of their sample. However, any dust-enshrouned RSGs will be so 
red that mimics the colors of AGBs. There is no effective way to distinguish very dusty 
RSGs from dusty AGBs with IR photometry solely, but they should have different features 



from the less-dusty counterparts ( iBoyer et al. 



201ll ). and they may also have different 



characteristics in the optical bands. 

In FigJ31 five outliers are marked by their ID numbers shown as Nos.3, 22, 27, 52 and 
86. Among them, No. 3, No. 52 and No. 86 are too faint to satisfy the lower limit of the 
Ks band luminosity. No. 22 exceeds the red borderline at J — Ks = 1.6 mag that may be 
caused by the heavy dust envelopes, and will be further investigated in MIR CMD since 
targets will have high luminosity in MIR if a large amount of dust exists. Although No. 27 
is a bit above the luminosity upper limit of a 30 M star, it is not excluded considering 
the uncertainty in the determination of the luminosity at the high end of the mass and the 
proximity to the borderline. 



Similarly, 



Boyer et al. 



(120111 ) divided the NIR CMD into several sub-regions for various 



types of ev olved stars as show n in Fig. HI where the J — Kg color cuts shown as solid 



lines follow 



Cioni et al. 



(12006a 



b|), with C-rich AGBs (C-AGBs) defined by K s < K0 and 



K s > K2, O-rich AGBs (O-AGBs) defined by K s < K0 and Kl < K s < K2, extreme-AGBs 
(x-AGBs) mostly by K s < K s -band ti p of the red giant branch (TRGB) (sho wn as dotted 



l ine) a nd J — Ks > 2 (see Figure 5 in 



Boyer et al 



201 ll ). Concerning RSGs, 



Boyer et al 



(120111 ) classified RSGs by restricting the the branch width to A(J — Ks) = 0.2 mag and 



leaving a 0.05 mag gap from the O-AG Bs shown as dashed lines in our diagram. The various 



regions defined by 



Boyer et al. 



(120 111 ) in the CMD are contoured in red, blue and purple 



solid lines for RSGs, O-AGBs and C-AGBs respectively. For a better view of distinguishing 
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the RSGs from AG Bs, all the AGBs, i.e. O-AGBs, C-AGBs and x-AGBs from Table 4 of 



Boyer et al. 



(1201 if ), are overplotted in this diagram as well. It is clear that the AGB stars 
are less luminous and redder than the RSGs. With these two restrictions, most of our 
targets apparently locate in the RSGs area (the branch width probably needs to be extended 
blueward) and only some dusty ones o r possible AGBs left out of the region. Almost all 



AGBs in Figure 5 of 



Boyer et al. 



(120111 ) are fainter than the limit of Ks = 9.5 mag (shown as 
dotted line) except one x-AGB. This confirms that our determination of the separative limit 
by converting from the bolometric magnitude is acceptable. The discrepancy between our 
and Boyer's results is that the magnitude range of their RSGs sample is twice of ours and 
has a downward tendency blueward to the Ks-TRGB. Meanwhile, our sample has a higher 
luminosity tip which indicate that more luminous sources are included. But according to 
the limit of Mboi = —7.1 mag and simple theoretical calculation of RSGs mass- luminosity 
relation (shown as dashed line in the diagram), Boyer's sample may contain lots of AGBs. 
In a word, the sample of Boyer has too low luminosity and too narrow and blue color index. 
In order to fit our sample, their region of RSGs should be moved redward 0.1 mag and 
upward 1.5 mag shown as a red dotted line in Fig. HI After such adjustment, the region 
looks reasonable for RSGs. The outliers in previous CMD are mostly in the AGBs area 
with No.3, N0.52 and No.86 possible being O-AGBs and No.22 belong an x-ABGs. 



2.1.2. Mid-infrared Color- Magnitude Diagrams and Two-Color Diagrams 



Most RSG candidates are already confirmed in the NIR CMDs. Considering RSGs are 
in the evolved phase and have large MLRs, MIR CMDs will shed light on the circumstellar 
dust of RSGs. S everal papers have pointed out that RSGs have circumstellar dust features 



from 8 to 12ixm (IHagen 



2008 



Verhoelst et al 



1978 



Skinner fc Whitmore 



1988 



Josselin et al. 



2000 



Ohnaka et al. 



20091 ). with the 9.7/im silicate feature prominent. 



13 



The color-magnitude diagram of [3.6]/J — [3.6] for all targets is shown in Fig. [5j The 



3.6 /iiii band is the most sensitive IRAC band comparable to the 2 MASS Ks band 



band is relatively free of dust and molecular emission/absorption (jVerhoelst et al. 



'he [3.61 



20091 ) 



The diagram of [3.6] versus J — [3 .6] reveals main 



where the region of AGBs set by 



Bonanos et al. 



y the property of the stellar photosphere, 



(120101 ) is shown as black solid line. The 



three outliers in previous NIR CMDs, Nos.3, 52 and 86, are within th e region of AGB st ars 



Boyer et al 



( 1201 If ) is 



so they may be AGB stars instead of RSGs. The region of RSGs by 
shown by red solid line. Also shown are the limits of the [3.6]-band magnitude derived from 
the mass of RSGs and the bolometric magnitude as dashed horizontal lines. No. 27 is still 
above the 3 line. The r ed ed ge of color index J — [3.6] is 2.6 mag, i.e. the reddest edge 



of AGBs by 



Bonanos et al. 



(120 101 ) since the d ust-producing rate s of RSGs are similar to 



AGB stars although RSGs are less numerous ( jBoyer et al. 



20111 ). Again, No. 22 is beyond 



this limit, indicating a really cold dust co mponent . Th e blue limit is J — [3.6] = 0.6 mag, 



i.e. the bluest color of RSGs region set by 



Boyer et al 



(l201l| ). With both the brightness 



and color limits, the outliers are the same as previous selection in the J/J — K s diagram, 
and it is also the same that most of our targets are settled in the RSGs region. 

Fig. [6] and Fig. [7] are CMDs of [8.0] versus K s — [8.0] and J — [8.0]. Following previous 
procedures, the [8.0] magnitude limits are set between [8.0] = 6.8 mag (for a mass of 30 
M Q ) and [8.0] = 9.9 mag (13 M Q ) in both diagrams, while these limits cannot be strict 
since the the [8.0] radiation comes not only from the photosphere but also the warm 
dust around the stars. The effect of dust is seen in the [8.0]/Ks — [8.0] diagram, where 
Ks — [8.0] becomes large with the increased [8.0] brightness, which can be fitted linearly by 
[8.0] = (-2.40 ± 0.09) (K s - [8.0]) + (9.59 ± 0.08) with a = 0.56 shown as a dotted line. 
Consistently, the red limit is set at Ks — [8.0] = 1.2 mag, the cross-point of the fitted line 
with the upper limit of luminosity at [8.0]. The blue limit is Ks — [8.0] = 0.1 mag since no 
targets bluer than this. It's interesting to notice that the 3645 massive stars in SMC cluster 
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in our two color limits (0.1 and 1.2 mag in Ks — [8.0]) respectively, while RSGs distribute 
quite evenly between the limits. As expected, there are more targets lying out of these 
restrictions, and most of the outliers are beyond the red limit, which can be explained by 
the dust emission. Besides the outliers Nos.3, 22, 27, 52, 86 already showing up in previous 
CMDs, there are new outliers Nos.2, 4, 23, 35, 55, 60, 77, 130, 14 all with K R - [8 ] > 1.2 
mag. In the [8.0]/ J — [8.0] diagram, the stellar classifications of iBoyer et all J2OII ) are 
shown as red solid lines. The brightness and color limits are set in the same way as in 
Fig. El except the linear fit is changed to be [8.0] = (-2.1 ± 0.07)(J - [8.0]) + (11.54 ± 0.10) 
with a = 0.57 that results in the limits of J — [8.0], 2.3 and 1.0 mag in the red and blue 



respectively. The das hed line at 



band (see Figure 4 in 



Boyer et al. 



8.0] = 8.4 mag is the upper limit of AGBs in the [8.0] 



20111) • The outliers, Nos.2, 3, 4, 22, 23, 27, 35, 52, 55, 77, 



86, 130, 140, are the same as before, No. 60 is no longer outlying and the new outlier No. 97 
is marked. 



T he MIPS [2 



band ( iBlum et al. 



band reflects the status of colder dust in compar ison with the [8.0] 



2006 



Bonanos et al. 



2009 



2010 



Bover et al.ll201lh . Fig. E shows the 



[24]/[8.0] — [24] CMD. The correlation of the brightness in [24] with the color index 
[8.0] — [24] is very tight since both reflect the radiation of the dust. The linear fit for this 
relationship is [24] = (-1.96 ± 0.05)([8.0] - [24]) + (9.03 ± 0.06) with a = 0.45. With 
the brightness limit in the [24] band between 5.0 and 10.3 mag for the mass range of 
RSGs, the color index [8.0] — [24] is correspondingly set to -0.4 and 2.1 mag. Considering 
the uncertainty in converting the mass to the brightness and the inclusion of the dust 
emission in the [24] band, the upper limit of the [24] band cannot be taken too strictly. 
The sources Nos. 27 and 131 are thus kept in the sample even though they lie above the 
upper limit line. In fact, the color index can also be influenced by the dust in addition to 
the photosphere, five targets, Nos. 22, 23, 35, 77, 130, with [8.0] — [24] redward of the limit, 
are neither excluded. On the contrary, their [8.0] — [24] indexes explain their deviation from 
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the majority in the NIR CMDs. For the same reason, Nos.3 and 52 previously classified as 
outliers are inside the limit box in this diagram. 

Similar to the CMDs, TCDs involving IR colors are good indicators to classify stellar 
populations. Combined CMD with TCD, the characteristics of RSGs in IR bands should 
be clear and distinguishable from other populations. 



Fig. [9] and Fig. [TT] are TCDs used in iBonanos et al.l (120101 ) and Fig. [10] is used in Paper 
I. The black solid and dashed lines represent the Bla ck Body (BB) a t vario us temperatures 
and a model wit h 3500 K BB plus 2 50 K dust from 



fsee Section 4 of 



Bonanos et ah 



Bonanos et al. 



(120101 ) respectively 



20091 for detailed description). Caution should be called 
that the spectra of RSGs dominated by emission lines make the blackbody not a perfect 
approximation to the atmosphere spectral energy distributions (SEDs). In Fig. [HJ most 
targets are located close to the point of the 3500 K BB. But there is still a discrepancy 
between the model and the targets that the targets are a bit bluer than the BB line, 
which may be due to the dust emission in the [8.0] band. The outliers are the same as in 
[8.0] /Kg — [8.0] and not marked in this diagram. 



Fig. [TOl shows the [3.6 
band by the CO lines (see 



[4.5] vs. [4.5] - [8 0] TCD. Due to the depression at the [4.5] 



Verhoelst et al. 



20091 ). most targets deviate from the BB model 



significantly. The obvious outliers in this diagram are Nos.4, 22, 23, 35, 85, 130. In Fig. [UJ 
the [3.6] - [8.0]/[8.0] - [24] TCD, many targets are below the 3500 K BB plus 250 K dust 
model line, which can be explained by the emission of silicate dust in the [8.0] band. 



2.1.3. Comparison of RSGs m SMC and LMC 



Paper I studied the RSGs in LMC. As the nearest two galaxies, both SMC and 
LMC present a possibility to spot the objects individually. With different metallicity, 
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Zsmc ~0.004 and Z 



LMC 



reveal the effect of metallicity on RSGs. 



0.008 (IGlatt et al.ll2010f). the com parison between them will 



Origlia et al 



( 119991 ) showed that higher metallicity 



brings about redder J — Kg by using instantaneous burst and continuous star formation 
models for the Galaxy and LMC (see Fig. 3 of their paper). The histogram of J — K s in 
Fig. [T2l confirms this conclusion. The peak value of J — Kg for RSGs in SMC is about 0.1 
mag bluer than in LMC. On the other hand, the bolometric magnitude hence the K§ band 
absolute magnitude is hardly de pendent on the metal licity with a difference between LMC 



and SMC smaller tha n 0.1 mag (jLevesque et al 



for LMC flAlves 



2006). With a distance modulus of 18.50 



20041 ) and 18.91 for SMC, the observed Ks magnitude is converted to the 
absolute magnitude Mk s in Fig. [131 where No. 22 with J — Ks = 2.79 mag is not included 
for an efficient view, and a shift of 0.1 mag towards red is added for the RSGs in SMC 
to compensate for the metallicity difference. It can be seen that most RSGs are mixed 
together in this diagram irrespective of their location in LMC or SMC. But comparing with 
33 targets in LMC, there are only 9 targets in SMC with Mk s < —10.84 mag corresponding 
to a mass of 25 M , which implies fewer high-mass RSGs in SMC. In addition, there are 
only 25 SMC targets with J — Ks > 1.2 mag while 48 such targets in LMC after the 0.1 
mag shift is corrected. This result confirm s that the dust produced by RSGs are less in the 



SMC than in the LMC (IBover et al 



2011) 



Given all these characteristics in color and brightness, the situation is complicate since 
some targets lack data in specific bands and some targets are within the limits in this 
criterion while beyond the limits in that criterion. In Fig. [HJ most targets that appear 
to be outliers in other diagrams actually are not far from the BB plus dust model, while 
Nos.4 and 130 are and thus excluded. Also excluded are Nos.3, 52 and 86 for their too low 
luminosities and no much dust to account for the low luminosity. In addition, Nos.34, 61 
and 132 lack almost all bands data. At last, 132 targets are identified as RSGs and kept for 
further analysis. 
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From the CMDs and TCDs, it can be seen that almost all massive sources from 



Bonanos et al. 



(120 lOf ) that locate in the same regions as RSGs are included in our sample. 
This indicates that our sample is nearly complete for the high-mass RSGs. But according to 
previous discussion, Boyer's sample may contain more low-mass RSGs, which needs further 



investigation since the lower mass limit of RSGs is somewhat uncertain. 



2.2. Period Determination 



The photometric data are collected from the ASAS and MACHO projects in the same 
way as in Paper I. Most of the photometry data comes from ASAS because the RSGs, 
except those heavily dust-enshrouded with extra visual extinctio n, easily get satura ted in 



2005|). The 



the MACHO observation due to their great intrinsic luminosity (IMassey et al. 
saturation limit is about 13 mag in the Kron-Cousins V band and about 14 mag in the 
R band (Kem Cook, private communicat ion) for the MACH O project. The photometric 



precision of MAC HO is about 0.02 mag (lAlcock et al. 



( iPojmanski 



2002|). 



19991 ) and of ASAS about 0.05 mag 



The way to handle the photometric data also follows Paper I. For ASAS, only the 
10-year long data in the standard Johnson V band is analyzed while the dropped I band 
data covers only about 500 days. The search radius is fixed to be 30" to retrieve the 
photometric data from the project web site and follow-up visual inspection of DSS image 
is performed to exclude confusions. Afterwards only the counterpart within 10" of the 
position is accepted (Grzegorz Pojmanski, private communication). The criteria exclude 22 
targets with poor measurement, 4 targets with more than two sources within the 10" circle 
and 1 target with no counterpart. With the exclusion of 27 sources from just identified 
132 candidates, 105 targets are left with a reliable association with the ASAS data. The 
original ASAS data is further processed for a better quality. First, a 3-sigma clipping is 
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applied to remove the deviating measurements (IHoaglin et al.lll983l ). Then the 



(Savitzky-Golay) polynomial smoothing filtering is performed to the light curve (IPress et al. 



east-square 



19921 ). Finally a 10-day bin is done to smooth out the short-term light variation. The 
process of handling the ASAS photometric data is shown in Fig. [TUfor No. 12 with original 
light curve, smoothed light curve and binned light curve respectively. 

For the MACHO data, the non-standard two-color photometric system are transformed 
to the standard Kron-Cousins V and R system in the same way as in Paper I. With a 
radius of 3" to search for the MACHO counterpart, 8 targets from the 132 candidates are 
found. One of them has the ASAS counterpart, three does not have and the others have 
poor ASAS photometry. Because of the low temperature and thick envelope, RSGs are 
usually much brighter in the R band than in the V band that makes R band easily get 
saturated and useless. In the process to convert the template magnitude into the standard 
system, the measurements with air mass larger than 2.0 are deleted. The measurements 
with photometric error bigger than 0.2 mag are also dropped. One more important thing 
concerning the MACHO data is the morphology of light c urve. Since the definition of R SGs 



2006|), the 



variables is semi-regular with amplitudes of about 1 mag (IKiss. Szabo fc Bedding! 
targets with regular periodic variation or/and amplitude larger than 1 mag are excluded. 
The left column of Fig. [15] exhibits the light curve of the three outliers in the MACHO 
database, i.e. Nos. 3, 52 and 86, with regular variation and large amplitude, which may be 
AGB stars with luminosity lower than RSGs shown in the CMDs. But there are 6 more 
sources in the 132 RSG candidates, Nos. 15, 22, 23, 60, 66 and 87, with similar light curve 
of regular variation or/and large amplitude, three of which are shown in the right column 
of Fig. [15j Excluding these 6 sources, only Nos. 46 and 59 are left as RSGs with reliable 
MACHO photometric data. Moreover, the sample of RSGs are reduced to 126 targets after 
excluding the 6 sources with regular or large amplitude light curves. 
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Combining the ASAS and MACHO data as a whole, there are 107 RSGs with available 
long-term visual photometry to derive the periods. The resource for the target is labeled by 
'M' for MACHO and 'A' for ASAS in Table [U 

RSGs are found to be variable for long time. But their variation is not regular that 
makes the period determination a tough task. On the other hand, a reliable determination 
of the period is the key to the P-L relation. In order to obtain the period as accurate as 
possible, four methods are tried to find the most consistent period. 

The first met hod used is the Phase Dispersion Minimization (PDM) method 



(IStellingwerl 



19781 ). while PDM2 takes place of the original version used in Paper I. PDM 
has many new features such as improved curve fits, suppressed subharmonics, and beta 
function statistics. This method is very useful in particular for data sets with gaps and 
non-sinusoidal variations. Since the light variation of RSGs is mostly non-sinusoidal and 
ASAS data have lots of gaps, it seems to be an appropriate technique. A sample of the 
PDM2 processing is shown in Fig. [16] for No. 12 (for comparison, the same target is taken 
as the example in discussing other methods). Two significant periods are present in the 
diagram with 1485 d (LSP) and 351 d respectively. The significance levels of 0.05 and 0.01 
are shown in dashed lines in the diagram. 



The second method to derive the period of light variation is Period04 (ILenz 



200J). 



Because the Long Secondary Period involved in the light variation of RSGs often exc eeds 



the t ime span of the observation and leads to an uncertain period by PDM, Period04 ( ILenz 



20041 ). based on the Fourier transform and least-square fitting, can make cross-confirmation. 
The Period04 processing of No. 12 is shown in Fig. [17] (see our Paper I for a detailed 
description of PDM and Period04 procedures). The results are similar to PDM2, with two 



2 http:/ /www. stellingwerf.com/ 
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periods at 361 d and 1553 d. 

The third method is newly taken into use this time, i.e the CLEAN method. This 
method deconvolves the spectral window function from the "dirty" Discrete Fourier 
ransform (DF T) spectrum by using a 1-D version of the iterative CLEAN algorithm 



( jHogbom 



1974]). During each iteration, the window function is centered at the amplitude 
peak. Then the full amplitude of this peak is inserted into the complex array of "CLEAN 
components" . To form a residual Fourier spectrum, the window function is scaled by the 
"CLEAN components" and subtracted from the input DFT. The process is repeated, with 
the residual spectrum from previous iteration used as the input spectrum. After several 
iterations, the "CLEAN components" is convolved with the Gaussian "beam", truncated 
at 5 sigma. Finally, the residual Fourier transform is added to this convolution to produce 
the CLEANed DFT spectrum. Through these steps, the information contained in the peak 
will be remained with the physical frequencies while their aliases or pseudo-aliases are 
decoupled. The example is shown in Fig. [18] The periods derived by the CLEAN method 
are 371 d and 1625 d. 

It can be seen that the above three methods derive quite consistent results, all with two 
periods and agreeing with each other within about 10%. But, there is still some uncertainty 
in the values, which may be c aused by the variation of the period itself. The Weighted 



Wavelet Z-Transform (WWZ) ((Foster 



19961 ) is used to evaluate the stability of the detected 
periods. The WWZ uses a standard Fourier fitting algorithm but with a sliding Gaussian 
weighting window whose width is determined by frequency and user-defined constant. 
The algorithm performs a sinusoidal wavelet to fit the data. The weight of each point is 
determined by the sliding window function, which means that the point close to the center 
of the window has heavy weights in the fitting. The window slides along the data set to 
represent the signal content of data as a function of time. The parameter c which weighs 
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the time resolution with frequency resolution is set to 0.005. The frequency range tested is 
from 0.0003 to 0.01 c/d with a step of 0.0001. An example is shown in Fig. [[Hi where the 
long period actually has a variation from about 1400 d to about 1600 d that explains the 
consistency of the PDM2, Period04 and CLEAN methods. 



3. Period-Luminosity Relation 



Based on the analysis of the time-series photometric data, the 126 targets which have 
been identified by CMDs, TCDs and light curve morphology are divided into five categories. 
(1) The first category includes 20 RSGs. They are either too bright and saturated in 
photometry or in a crowded stellar field with a close companion unresolvable. These stars 
are not taken into account in further investigation. (2) The second category includes 55 
RSGs. They are irregular variables with so complex light curve that no reliable period can 
be determined. They are neither considered for further study of the P-L relation. The 
other 51 RSGs are semi-regular or LSP RSGs which are involved in the determination of 
the P-L relation. (3) The third category includes 6 RSGs, being semi-regu lar variables with 



period statistical significance less than or equal to 0.05. They are all from 



Massey Olsen 



(120031 ). (4) The fourt h category includes 15 R SGs with 



period, 1 targe t from 



Levesque et al 



Catchpole &: Feast 



(120061 ) and 1 fro mMassey et al 



1981), 11 fr om 



;SP and distinguishab 



Massev & Olsen! (I2003f ). 2 from 



e short 



(120071 ). (5) The fifth category includes 30 



RSGs with LSP but no distinguishable short period. The relatively high-accuracy period 
can only be achieved for the third and fourth categories, i.e. the 21 RSGs with semi-regular 
variation or with distinguishable short period variation. It's them that are taken into 
account when deriving the P-L relation afterwards. Their periods are listed in Table El with 
also the period averaged over three methods and the dispersion most of which are shorter 
than 10 days. For an intuitive view, examples of irregular, semi-regular and LSP RSGs 
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variation as well as the two light curves from the MACHO database are shown in Fig. 

We also calculate the relation of period and amplitude for the RSGs in Table [3J Since 
there are two groups of period and amplitude derived by PDM and Period04 respectively, 
they are combined to make a robust linear fitting. The result is 

AV = (0.76 ± 0.14) x log P - (1.83 ± 0.24) 
with a = 0.10, where AV is the amplitude in the V band and P is the period. Except few 
outlying points, most points roughly obey the relation. This is consistent with the general 
tendency of variables that the longer the period the greater the amplitude. The distribution 
of the amplitude and the period is shown in Fig. 



With the period of variation determined, luminosity is the other key parameter in 
deriving the P-L relation and deserves as important attention as period. Historically, V 
band is the most common observed and used as the indicator of luminosity, but visual band 
suffers interstellar and/or circumstellar extinction severely. 

Fortunately, the extinction influence is much smaller in near- or mid-infrared bands, 
e.g. the extinctio n in the GLIPMSE/IRAC bands is only about 25th of the visual extinction 



( iGao et al. 



20091 ). Therefore, the infrared bands, 2MASS JHK S , IRAC and MIPS 24 
[xm band, are taken as the luminosity indicator while the V band is included only for 
comparison. 



3.1. K s -band Period-Luminosity Relation 

In Paper I, the P-L relation is determined for the LMC RSGs based on 47 semi-regular 
and LSP RSGs with distinguishable short period. Here a similar work is extended to the 
SMC, based on 6 semi-regular and 15 LSP RSGs with distinguishable short period listed in 
Table [31 Fig. [22] shows the robust linear fitting for the Ks-band P-L relation of SMC RSGs, 
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where filled and unfilled squares denote the semi-regular and LSP RSGs respectively. The 
fitting takes a form of m^ = a log P + b and the result is 

K s = (-3.28 ± 0.39) x logP + (17.25 ± 0.64) 
With a dispersion of o = 0.27, the relation is very tight. 

The P-L relation h as be en calculated for RSGs in the Galaxy and M33 by 



Pierce. Jurcevic fc Crabtree! (120001 ) (52 



Kiss. Szabo &: Bedding! (120061 ) (13 targets) and 
targets) respectively. For independence and consistency, we re-calculate the P-L relation 
in the Kg bands with some modification because the photometric data in infrared is only 
available in the Kg band for RS Gs in the Galaxy and M33. Th e absolute Kg-band magnitude 



of Galactic RSGs is taken from 



Kiss. Szabo fc Bedding! (120061 ) (see section 4.1 in their paper 



for detail description), but with ad dition of 3 targets, GP Cas, YZ Per and V500 Cas from 



„ ^^-r.v.-w^,, . 

Pierce. Jurcevic & Crabtree 


(2000) 


2004 


) since they are all members 


Pierce. Jurcevic & Crabtree 


(2000) 



Kinman. Mould fc Wood! (119871 ) 



who investigated red long-period variables including both RSGs and AGBs. For the sake 
of simplicity, only the lower Kg band magnitude cut-off is set as Kg = 14.86 corresponding 
to the lower mass limit of RSGs (9 M ) to the original dataset of 



Kinman. Mould fc Wood 



(119871 ) for 54 targets. This excludes 14 targets and keeps 40. The distance moduli of LMC, 



SMC and M33 used to convert from t 



18.50, 18.91 and 24.93 dBonanos et al 



re ob s erved Kg to t he absolute magnitude M Ks are 



2006 



U et al. 



20091 ) respectively. 



After these modifications, a universal linear fitting of the Kg band P-L relation is 
performed for RSGs in the Galaxy, LMC, SMC and M33. The results are shown in Fig. 
and Table. HI It can be seen that a tight P-L relation exists for RSGs in all these galaxies. 
In particular, the relation is quite consistent for the three extra-galaxies LMC, SMC and 
M33, while the Galaxy has a slight different slope and zero point which may be caused by 
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the distance error. Combining all the data together, the linear fit results in 



M 



-3.68 ± 0.18) x logP + (-0.62 ± 0.30) 



with a = 0.31. Also show n in Fig. [2~4l are the P-L rela t ions o f RSGs from 



Catchpole k Feast 



(2000) and 



(119801 ) and 



(Il98lh. 



Kinman Mould k Wood! (119871 ) 



Kiss. Szabo k Bedd 



Catchpole k Feast 



Feast et al. 



(Il980h . 



Pierce. Jurcevic k Crabtree 



nd (20061 ) by different line styles. The results of 



Feast et al. 



(119811 ) obviously deviate from ot hers, which may be caused 



by the i r impure sample contaminated by A GBs . The results of 



(1987) 



Pierce. Jurcevic fe Crabtred (120001 ) and 



Kinman. Mould k Wood 



Kiss. Szabo k Bedding! (120061 ) have quite 



close slope and intercept as ours. 



Within the SMC , 



and 



Soszyhski et al. 



Kiss k Beddingl (120031 ) 



Ita et al. 



(J2004J) 



Soszyhski et al. 



(120071 ) 



(120111 ) analyzed the P-L relations of LPVs found by the Optical 



Gravitational Lensing Experiment (OGLE) survey. The comparison shows that our result 



matches well with the extension of the AG B a? sequence o 



Red Giants (OSARGs) variables defined by 



Soszyhski et al. 



;he O GLE Small Amplitude 



(120071 ). For the OSARGs, the 



sequence 02 is the second sequen ce (from longest to shortest periods) that is similar to the 



sequence B of 



giants. But 



Wood et al. 



Kiss k Beddingl (120031 ) and 



(|1999l) a nd i dentified a s the 



Ita et al. 



irst overtone radial pulsating red 



(120041 ) noticed an additional sequence 



(labeled C 7 ) between Wood's sequences B and C. If sequence C is identified as the 
fundamental mode of pulsation, sequence C would be the first overtone and sequence B 
would be the second overtone, which makes it suspicious to identify the pulsation mode of 
RSG as the first overtone pulsation from the P-L relation. With the help of the theoretical 
model prediction, the pulsation mode of RSGs will be clearer as will be discussed in Section 
3.3. It is also noticed that there are only four targets with the period shorter than 250 d, 
which may be caused by two reasons. One is the sample selection effect that the lower-mass 
RSGs with the shorter periods are missed. This can be seen also in CMDs of this paper 
that most targets have theoretical mass larger than 13 M Q . The other is the ambiguity of 
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the relation between pulsation and convection for RSGs of different masses. Usually, the 
longer the period the greater the amplitude, vice versa. However, due to the existence of 
large convection cells in the surface of RSGs, they could produce unpredictable variation to 
even "overwhelm" the pulsation if the period is short enough, viz amplitude small enough. 
Thus no periodic variation can be observed. It can be seen in Fig. [21] of this paper and 
Fig. 10 in Paper I that the amplitude decreases to an undetectable level compared to the 
observational error if period shorter than about 250 d (logP < 2.4). 

3.2. Other infrared P-L relations 

The photometric data are available in other infrared bands than the Ks band by space 
telescopes such as Spitzer. The study of RSGs in LMC already found that the P-L relations 
exist in other infrared bands. The change of the brightness with period is thus investigated 
in the infrared bands, namely the 2MASS JHKs bands, the Spitzer /IRAC bands and the 
Spitzer /MIPS 24 /zm band, as well as the visual band V. Similar conclusions are drawn for 
RSGs in SMC as in LMC. The V band magnitude does not have clear correlation with the 
period because of the severe effect of extinction, in particular the RSGs in SMC are mostly 
located in the star-forming region as can be seen in Fig. [TJ No effort is put into fitting the 
P-L relation in visual. On the contrary, there is little effect of extinction on the brightness 
in the infrared bands. Moreover, the amplitude of variation in infrared is much smaller than 
in visual. Thus, the one-epoch photometric result can be a good and stable indicator of the 
luminosity of the object. One just needs to be careful that the infrared brightness partly 
comes from the circumstellar dust instead of the stellar photosphere, which becomes more 
significant at the longer wavelengths, i.e. the [8.0] and [24] bands in our study. Nevertheless, 
a linear fitting is performed in all the related infrared bands in the form of M^ = a log P + b 
for LMC and SMC separately and for their combination. The slope a and intercept b with 
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dispersions are listed in Table [51 and the fitting diagrams are shown in Fig. [25] and Fig. 



The diagrams tell that the P-L relationship does exist in the RSGs in SMC, the same 
as in LMC, in all infrared bands with the dispersion depending on the band. In the V band, 
the dispersion is so large that no relation is intended to determine. But the LMC RSGs are 
on average fainter than the SMC RSGs in the V band, which can be caused by the higher 
metallicity in LMC for more circumstellar dust extinction as can also be seen in the [8.0]- 
and [24]-related diagrams. In the J, [8.0] and [24] bands, the dispersions are the moderate, 
and the P-L relations are reasonably good. In the H, Ks, [3.6], [4.5] and [5.8] bands, the 
dispersions are very small, and the P-L relations are very tight and reliable. The best 
relationship occurs in the [4.5] band, the same as for LMC in Paper I. This is reasonable 
since the radiation mainly comes from the stellar photosphere if the targets are cold and 
have no thick dust envelope and the extinction is very small in these bands. So the IRAC 
bands as well as the Ks band are recommended when adopting one-epoch photometry as 
the luminosity indicator. 

Long secondary period is another period for those RSG with LSP. All the LSP targets 
have relative long period and small amplitude compared to the semi-regular RSGs, which 
means a different origin. Many models such as radial and non-radial pulsation, binary, 



stellar spot mode 



flWood fc Nicholls 



s, have been su 



2009 



ggeste d but none of them receive s general acceptance 



Nieetal 



20101 ). Recently, 



Stotherd ( 120101 ) has demonstrated that 



giant convection cells in the envelopes of massive red supergiants turn over in a time 
comparable in order of magnitude to the observe d long secondary periods in these stars 



Stothers fc Leung! (119711 ). He has applied 



according to a theory proposed some years ago by 
the theory successfully to the well observed RSGs of Betelgeuse, Antares and M-type giant, 
while more samples may need to prove this model. It will be interesting to know if LSP 



is related to luminosity. The amplitude and period of 45 SMC LSP are calculated in the 
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same way as for the semi-regular RSGs by three methods. Because the observation spans 
only about 3000 days, 14 targets with period longer than this are discarded, while the 
results for the left 31 targets are listed in Table O Fig. |2"T1 shows their absolute magnitude 
in the Ks band versus the period, together with the LSP of RSGs in LMC and the short 
periods of RSGs in the four galaxies dis cussed earlier. S u perpo sed in this diagram are the 



P-L relations of LPVs (solid 



(points) from 



Soszyhski et al 



i nes) 



rom ISoszyhski et al. 



( 120071 ) and various LPVs targets 



(1201 if ) in the SMC, where the sym bols are the same as th e 



original papers. With the addition of more a nd fainter objec t s in 



there are small shifts of the P-L relations by 



Soszyhski et al. 



Soszyhski et al 



(2011). 



( 120071 ) from the locations of 



the objects. There is no clear P-L relation for these long secondary period. The major 
obstacle of LSP investigation is the need of very long-term observation and the irregularity 
of the variation. With the continuation of the ASAS project, a more reliable determination 
of the period and analysis of the mechanism are possible. 



3.3. Pulsation mode 



The exact mode of pulsation of RSGs is unclear although the radial modes are 
commonly accepted. In order to understand the pulsation mode of RSGs, the P-L 
relation derived from the observation is compared with the theoretical P-L relation. The 
observed bolometric magnitude needs to be converted to stellar luminosity to conveniently 
compare with the mo d el. The reliable absolu te b olometric magni t udes are collected from 



Levesque et al. 



(120061), 



Levesque et al. 



new RSGs effective temperature scales (ILevesque et al 



(I2007h and 



Levesque et al. 



(120091 ) which applied the 



20051 ) to LMC and SMC to derive 



M bol that corrected the defect of the old effective temperature. The adopted M bo i is the 
value derived from spectral fitting other than (V— K) which can be significantly affected 
by the extinction in the V band. This selection of Mboi picks out only 13 targets in Table. [3] 
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in SMC and 13 targets in LMC. 



The t heoretical P-L re lations of RSGs were predicted for the first three radial pulsation 
modes by iGuo fc Lil (120021 ) at three typical masses, 15, 20 and 30 M Q based on a linear 



non-adiabatic model. In particular, they considered the rela tions under differen t metallicity 



that covers the cases of SMC, LMC and Galaxy. Similarly, 



Heger et al. 



(119971 ) computed 



the P-L relations of RSGs from linear non-adiabatic models for four masses, 1 0, 12, 15 



and 2 M , but without discrimination of the metallicity. The P-L relations by 



Guo&Li 



( 120021 ) are derived by simulating the evolutionary tracks of the stars in the H-R diagram 
with different chemical composition and a constant mixing- length parameter a = 1.0. Their 
results show that period increases almost monotonically in the RSG phase (see Table 3 
and Figure 4-12 in their paper). Our results are compared with their P-L diagrams for 
Z = 0.005, Z = 0.01 and Z = 0.02, as the best approximations to the SMC (Z = 0.004), LMC 
(Z = 0.008) and Galaxy (Z = 0.02) respectively. In Fig. [28] for the SMC case, our targets 
are mostly located in the range of mass between 20 and 30 M Q . From the observed fact that 
luminosity of RSG increases with period and the theoretically predicted P-L relations for 
three typical masses, our targets lie in the region of the theoretical first-overtone mode. The 
conclusion is similar in Fig. [29] for Z = 0.01 for the LMC RSGs, with a slight shift to the 
region of the fundamental mode. Two objects may even belong to the fundamental-mode 
region. The shift in Fig. [301 for the Galaxy case is more significant. Half of the RSGs in the 



Galaxy from 



Kiss. Szabo & Bedding! ( 120061 ) fits the P-L relations at the fundamental-mode 



a nd the other half at the first-overtone mode. When compared with the theoretical models 



of 



Heger et al. 



(119971 ) in Fig. [31] the conclusion is the same, i.e. most objects are consistent 



with the predicted P-L relation at the first-overtone mode except that two in LMC and half 
in the Galaxy are close to the fundamental-mode relation. Indeed, the identification of the 
first overtone mode is already implied in that the P-L relations of RSGs in both LMC and 
SMC agree very well with the extension of the AGB sequence corresponding to the first 



-29 - 



overtone radial mode of OSARGs (jSoszyriski et al. 



20071 ). However, as mentioned earlier 



and seen in Fig. [27J the P-L relation of RSGs are not perfectly matched with the AGB 
a 2 sequence with the newly accumulating data while lies a kind of between the first- and 
second-overtone modes. This inconsistency may indicate that RSGs have their own P-L 
sequence, which deserves further investigation. It may be too exaggerated to argue that the 
pulsation mode is related to the metallicity only from these three cases of SMC, LMC and 
the Galaxy, while this deserves further study. Another approach to clarify the pulsation 
mode relies on the accurate determination of the absolute bolometric magnitude for more 
complete sample of RSGs, especially those in the Galaxy. 



4. Summary 

The P-L relation of RSGs has aroused much interest since 1980s due to the evidence 
that such luminous and pulsatile objects can lead to an extension of stellar standard candle. 
Using long-term optical survey and large-scale infrared survey, we have the opportunity to 
calibrate the P-L relation in several infrared bands for nearby Local Group galaxies and 
identify the pulsation mode with the help of the theoretical calculation. 

The preliminary sample of RSGs in SMC is consisted of 140 stars assembled from 
previous works. The infrared counterparts are chosen in the 2MASS, Spitzer /IKAC or 
MIPS catalogs by a one-arcsec position consistency and identified by the brightness and 
colors in the CMDs and TCDs. These procedures eliminate 5 targets for inconsistent color 
or luminosity and 3 targets for lack of brightness information, leaving 132 targets as the 
RSG candidates. 

The time-series photometry data are taken from the databases by the ASAS and 
MACHO projects and further processed. The period of light variation is derived by the 
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PDM, Period04, CLEAN and WWZ methods. According to the morphology of the light 
curve, 6 more sources are further identified as AGBs and removed from the RSGs sample. 
Based on the properties of light variation, the 126 RSGs are classified into five categories: 
20 with poor photometry, 55 with no reliable period, 6 with semi- regular variation, 15 with 
LSP and distinguishable short period and 30 with only LSP. 

For the semi-regular variables and the LSP variables with distinguishable short period, 
the Ks band P-L relation is analyzed and determined. This relation is compared with those 
in the Galaxy, LMC, and M33. It is found that the RSGs in these galaxies obey similar P-L 
relation except the Galaxy. In addition, the P-L relations in the infrared bands, namely the 
2MASS JHK S and Spitzer/IRAC and Spitzer /MIPS 24 jim band, are derived with high 
reliability. The best P-L relation occurs in the Spitzer/IRAC [3.6] and [4.5] bands. No clear 
P-L relation is found for the LSP RSGs. 

From the P-L relation, the mode of pulsation of RSGs in SMC is discussed and found 
to be very possibly in the first overtone radial mode. When compared with the Galaxy and 
LMC, it seems there is a tendency that the pulsation mode depends on the metallicity. 

The authors thank Prof. L. L. Kiss for kindly providing the data of the Galactic RSGs. 
This work is supported by the China's NSFC projects Nos. 10973004 and 11173007. BWJ 
thanks the John Templeton Foundation for support. 
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Table I. Infrared Catalog of the 140 SMC RSGs candidates 



No. Name R.A.(J2000) Decl.(J2000) J H K s [3.6] [4.5] [5.8] [8.0] [24] Spectral Type a Reference 13 Source c 



1 


HV11223 


00 


32 


01 


.61 


-73 


22 


34 


.7 


11.166 


10.276 


9.971 


2 


HV1349 


00 


41 


21 


.50 


-72 


50 


15. 





11.202 


10.446 


9.946 


3 


HV1366 


00 


42 


49. 


.87 


-72 


55 


11. 


.4 


12.232 


11.431 


11.154 


4 


HV1375 


00 


42 


52 


.30 


-73 


50 


51. 


.0 


11.300 


10.381 


9.832 


5 


[M2002]SMC005092 


00 


45 


04 


.56 


-73 


05 


27. 


.4 


9.208 


8.377 


8.071 


6 


[M2002]SMC008324 


00 


47 


16 


.84 


-73 


08 


08. 


4 


10.456 


9.685 


9.507 


7 


[M2002]SMC008367 


00 


47 


18 


.11 


-73 


10 


39. 


3 


10.080 


9.261 


9.040 


8 


[M2002]SMC008930 


00 


47 


36 


.94 


-73 


04 


44. 


3 


9.466 


8.601 


8.319 


9 


[M2002]SMC009766 


00 


48 


01 


.22 


-73 


23 


37. 


.5 


10.208 


9.402 


9.213 


10 


[M2002]SMC010889 


00 


48 


27 


.02 


-73 


12 


12 


3 


8.778 


7.990 


7.774 


11 


[M2002]SMC011101 


00 


48 


31 


.92 


-73 


07 


44. 


.4 


10.316 


9.468 


9.237 


12 


[M2002]SMC011709 


00 


48 


46 


.32 


-73 


28 


20 


.7 


9.510 


8.770 


8.530 


13 


[M2002]SMC011939 


00 


48 


51 


.83 


-73 


22 


39 


3 


9.629 


8.864 


8.606 


14 


[M2002]SMC012322 


00 


49 


00 


.32 


-72 


59 


35. 


.7 


9.280 


8.456 


8.276 


15 


CV78 


00 


49 


03 


.84 


-73 


05 


19 


.9 








16 


[M2002]SMC012572 


00 


49 


05 


.25 


-73 


31 


07. 


.8 


9.204 


8.478 


8.258 


17 


[M2002]SMC012707 


00 


49 


08 


.23 


-73 


14 


15 


.5 


10.365 


9.560 


9.333 


18 


[M2002]SMC013472 


00 


49 


24 


.53 


-73 


18 


13. 


.5 


9.111 


8.485 


8.301 


19 


[M2002]SMC013740 


00 


49 


30 


.34 


-73 


26 


49 


.9 


10.478 


9.623 


9.387 


20 


[M2002JSMC013951 


00 


49 


34 


.42 


-73 


14 


09. 


9 


10.084 


9.262 


9.055 


21 


[M2002]SMC015510 


00 


50 


06 


.42 


-73 


28 


11. 


1 


9.590 


8.760 


8.535 


22 


IRAS F00483-7347 


00 


50 


07 


.19 


-73 


31 


25. 


.2 


11.431 


9.777 


8.639 


23 


IRAS X0048-731 


00 


50 


30 


.63 


-72 


51 


29. 


9 


10.338 


9.394 


8.783 


24 


[M2002]SMC017656 


00 


50 


47 


.22 


-72 


42 


57. 


2 


9.839 


9.080 


8.921 



9.749 


9.878 


9.683 


9.557 


9.205 





CF81 


A 


9.284 





8.874 


8.255 


6.230 





CF81 


A 


10.733 


10.594 


10.539 


10.283 


9.844 





CF81 


M 


8.715 


8.365 


7.835 


7.204 


5.388 





CF81 


— 


7.765 


7.733 


7.491 


7.081 


5.185 


M2 I 


L06 


A 


9.372 


9.462 


9.367 


9.267 


9.608 


K0: I 


M03 


A 


8.838 


8.867 


8.708 


8.501 





K7 I 


M03 


A 


8.039 


8.043 


7.814 


7.641 


7.235 


MO I 


M03 


A 


9.056 


9.219 


9.069 


8.961 


8.644 


K7 I 


M03 


A 


7.571 


7.682 


7.463 


7.285 


5.628 


M0 I 


M03 


A 


8.978 


9.091 


8.945 


8.879 


8.836 


K2.5 I 


M03 


A 


8.345 


8.471 


8.299 


8.205 


7.479 


K5-M0 I 


M03 


A 


8.296 


8.412 


8.116 


7.894 


5.936 


K2 I 


M03 


A 


8.045 


8.126 


7.937 


7.799 


6.522 


K5-M0 I 


M03 


A 


9.311 


9.193 


8.955 


8.760 


7.985 




WBF83 


M 


8.125 


8.236 


8.159 


8.118 


7.996 




M03 


A 


9.140 


9.204 


9.077 


8.954 


8.845 




M03 


A 


8.084 


8.013 


7.823 


7.484 


7.242 


K7: I 


M03 


A 


9.241 


9.363 


9.219 


9.145 




K3 I 


M03 


A 


8.870 


9.023 


8.859 


8.776 


8.485 


K3 I 


M03 


A 


8.280 


8.401 


8.197 


7.965 


6.701 


K5 I 


M03 


A 


6.620 


6.309 


5.507 


4.790 


2.591 




V08 


M 


8.643 


8.212 


8.030 


7.030 


4.740 




V08 


M 


8.776 


8.861 


8.737 


8.647 


8.737 


KO-5 I 


M03 


A 



Table 1 — Continued 



No. 


Name 


R.A.(J2000) 


Decl.(J2000) 


J 


H 


K S 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


Spectral Type a 


Reference 13 


Sourc 


25 


HV11295 


00 


50 


48.40 


-72 


52 


29.0 


10.849 


10.004 


9.580 


9.267 


9.120 


8.948 


8.804 


7.910 




CF81 


A 


26 


[M2002]SMC018136 


00 


50 


56.01 


-72 


15 


05.7 


8.904 


8.107 


7.849 


7.631 


7.819 


7.609 


7.452 


5.850 


MO I 


L06 


A 


27 


[M2002]SMC018592 


00 


51 


03.90 


-72 


43 


17.4 


8.385 


7.671 


7.425 


6.885 




6.706 


6.304 


4.605 


KO-2 I 


M03 


A 


28 


[M2002]SMC019551 


00 


51 


20.23 


-72 


49 


22.1 


10.362 


9.593 


9.382 


9.239 


9.347 


9.233 


9.152 


9.162 


K2 I 


M03 


A 


29 


[M2002]SMC019743 


00 


51 


23.28 


-72 


38 


43.8 


9.678 


8.899 


8.730 


8.451 


8.466 


8.221 


7.947 


6.372 


K5 I 


M03 


A 


30 


[M2002]SMC020133 


00 


51 


29.68 


-73 


10 


44.3 


9.220 


8.410 


8.155 


7.953 


7.957 


7.693 


7.383 


5.509 


MO I 


M03 


A 


31 


[M2002]SMC020612 


00 


51 


37.57 


-72 


25 


59.5 


10.451 


9.737 


9.569 


9.412 


9.511 


9.394 


9.302 


9.215 


K5 I 


M03 


A 


32 


[M2002JSMC021362 


00 


51 


50.25 


-72 


05 


57.2 


9.862 


9.017 


8.824 


8.674 


8.749 


8.573 


8.408 


8.210 


K5 I 


L06 


A 


33 


[M2002JSMC021381 


00 


51 


50.46 


-72 


11 


32.2 


9.967 


9.206 


8.999 


8.797 


8.844 


8.708 


8.609 


8.440 


K5 I 


L06 


A 


34 


HV11303 


00 


52 


09.00 


-71 


36 


22.0 




















CF81 


A 


35 


ISO-MCMS J005212.9- 730852 


00 


52 


12.82 


-73 


08 


52.8 


11.254 


10.468 


9.917 


8.800 


8.487 


8.268 


7.467 


4.734 




V10 


A 


36 


[M2002]SMC023401 


00 


52 


25.36 


-72 


25 


13.3 


10.358 


9.627 


9.425 


9.303 


9.370 


9.321 


9.192 


9.095 


Kl I 


L06 


A 


37 


[M2002]SMC023463 


00 


52 


26.51 


-72 


45 


15.6 


9.674 


8.956 


8.730 


8.328 


8.341 


8.101 


7.926 


7.643 


KO-5 I 


M03 


A 


38 


[M2002]SMC023700 


00 


52 


30.69 


-72 


26 


46.8 


10.424 


9.656 


9.440 


9.325 


9.477 


9.341 


9.267 


9.106 


KO-2 I 


M03 


A 


39 


[M2002]SMC023743 


00 


52 


31.49 


-72 


11 


37.3 


10.304 


9.570 


9.411 


9.230 


9.314 


9.197 


9.077 


9.100 


K2 I 


L06 


A 


40 


[M2002]SMC025550 


00 


53 


02.85 


-73 


07 


45.9 


10.500 


9.677 


9.453 


9.301 




9.303 




9.044 


K2 I 


M03 


A 


41 


[M2002]SMC025879 


00 


53 


08.87 


-72 


29 


38.6 


9.359 


8.624 


8.447 


8.139 


8.146 


8.021 


7.835 


6.995 


MO I 


M03 


A 


42 


[M2002]SMC025888 


00 


53 


09.04 


-73 


04 


03.6 


9.262 


8.454 


8.217 


7.983 


8.043 


7.905 


7.713 


6.743 


K5-7 I 


M03 


A 


43 


[M2002]SMC026402 


00 


53 


17.81 


-72 


46 


06.9 


10.276 


9.487 


9.271 


9.068 


9.178 


9.003 


8.856 


8.278 


KO-2 I 


M03 


A 


44 


[M2002]SMC026778 


00 


53 


24.56 


-73 


18 


31.6 


10.188 


9.433 


9.222 


9.037 


9.150 


8.992 


8.892 


8.631 


K2 I 


M03 


A 


45 


[M2002]SMC027443 


00 


53 


36.44 


-73 


01 


34.8 


9.634 


8.873 


8.640 


8.388 


8.521 


8.347 


8.221 


7.360 


K7 I 


M03 


A 


46 


HV11329 


00 


53 


39.40 


-72 


52 


39.0 


10.911 


9.987 


9.647 


9.493 


9.562 


9.416 


9.302 


9.076 




CF81 


M 


47 


[M2002]SMC027945 


00 


53 


45.74 


-72 


53 


38.5 


10.434 


9.820 


9.619 


9.496 


9.566 


9.443 


9.302 


9.214 


K3-5 I 


M03 


A 


48 


[M2002]SMC030135 


00 


54 


26.90 


-72 


52 


59.4 


10.335 


9.690 


9.492 


9.236 


9.408 


9.161 


8.892 


8.048 


K2 I 


M03 


A 



Table 1 — Continued 



No. 


Name 


R.A.(J2000) 


Decl.(J2000) 


J 


H 


K S 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


Spectral Type a 


Reference b 


Source' 


49 


[M2002JSMC030616 


00 


54 


35.90 


-72 


34 14.3 


9.308 


8.531 


8.341 


8.125 


8.199 


8.017 


7.768 


7.333 


K2 I 


M03 


A 


50 


[M2002]SMC032188 


00 


55 


03.71 


-73 


00 36.6 


9.728 


9.052 


8.861 


8.765 


8.922 


8.782 


8.662 


8.358 


K5 I 


M03 


A 


51 


HV12122 


00 


55 


18.01 


-72 


05 31.9 


10.600 


9.714 


9.361 


8.911 


8.834 


8.558 


8.419 


8.376 




WBF83 


A 


52 


HV1645 


00 


55 


19.76 


-73 


14 42.2 


12.663 


11.775 


11.421 


10.812 


10.674 


10.492 


10.290 


10.010 




CF81 


M 


53 


[M2002]SMC033610 


00 


55 


26.82 


-72 


35 56.2 


9.803 


9.014 


8.762 


8.558 


8.604 


8.409 


8.195 


6.445 


K7 I 


M03 


A 


54 


[M2002]SMC034158 


00 


55 


36.58 


-72 


36 23.6 


9.877 


9.123 


8.909 


8.727 


8.737 


8.543 


8.365 


7.185 


K2 I 


M03 


A 


55 


HV 838 


00 


55 


38.21 


-73 


11 41.1 


10.608 


9.860 


9.440 


8.369 




8.168 


8.148 


7.513 




CF81 


A 


56 


[M2002]SMC035231 


00 


55 


55.10 


-72 


40 30.4 


9.995 


9.482 


9.343 


9.167 


9.142 


9.000 


8.935 


8.782 


K2 I 


M03 


A 


57 


[M2002]SMC035445 


00 


55 


58.84 


-73 


20 41.4 


9.976 


9.177 


8.980 


8.834 


9.038 


8.855 


8.765 


8.673 


Kl I 


L06 


A 


58 


[M2002]SMC037994 


00 


56 


43.55 


-72 


30 15.0 


9.779 


8.946 


8.691 


8.518 


8.599 


8.451 


8.354 


8.154 


K7 I 


M03 


A 


59 


HV11366 


00 


56 


54.77 


-72 


14 08.6 


11.066 


10.154 


9.870 


9.678 


9.795 


9.608 


9.518 


9.458 




CF81 


M 


60 


HV1719 


00 


57 


14.48 


-73 


01 21.3 


10.780 


10.117 


9.798 


9.091 


8.979 


8.734 


8.609 


7.913 




CF81 


A M 


61 


[M2002]SMC041778 


00 


57 


56.45 


-72 


17 33.3 
















8.325 


K7 I 


M03 


A 


62 


[M2002]SMC042319 


00 


58 


06.61 


-72 


20 59.8 


10.178 


9.364 


9.137 


8.987 


9.112 


8.955 


8.873 


8.637 


K2-5 I 


M03 


A 


63 


[M2002]SMC042438 


00 


58 


08.71 


-72 


19 26.7 


10.383 


9.596 


9.362 


9.192 


9.321 


9.206 


9.152 


9.078 


K2 I 


M03 


A 


64 


[M2002]SMC043219 


00 


58 


23.30 


-72 


48 40.7 


10.173 


9.368 


9.114 


8.886 


9.063 


8.899 


8.838 


8.536 


K2 I 


M03 


A 


65 


[M2002]SMC043725 


00 


58 


33.21 


-72 


19 15.6 


10.595 


9.774 


9.558 


9.391 


9.542 


9.416 


9.338 


9.239 


K5 I 


M03 


A 


66 


HV12149 


00 


58 


50.17 


-72 


18 35.6 


9.963 


9.022 


8.609 


8.548 


8.381 


8.192 


7.889 


6.878 




CF81 


M 


67 


[M2002]SMC045378 


00 


59 


07.16 


-72 


13 08.6 


10.039 


9.205 


9.005 


8.807 


8.956 


8.814 


8.697 


8.694 


K3 I 


M03 


A 


68 


[M2002]SMC045850 


00 


59 


16.90 


-72 


25 10.9 


10.150 


9.366 


9.199 


9.004 


9.131 


8.996 


8.929 


8.807 


KO-5 I 


M03 


A 


69 


[M2002]SMC046497 


00 


59 


31.33 


-72 


15 46.4 


9.432 


8.589 


8.308 


8.160 


8.145 


7.986 


7.788 


6.212 


K5-M0 I 


M03 


A 


70 


[M2002]SMC046662 


00 


59 


35.04 


-72 


04 06.2 


9.377 


8.627 


8.345 


8.160 


8.209 


8.001 


7.799 


6.142 


MO I 


M03 


A 


71 


[M2002]SMC046910 


00 


59 


40.58 


-72 


20 55.9 


10.227 


9.539 


9.330 


9.181 


9.258 


9.135 


9.021 


8.943 


K3-5 I 


M03 


A 


72 


[M2002]SMC047757 


01 


00 


00.63 


-72 


19 40.2 


9.294 


8.502 


8.240 


8.084 


8.104 


7.894 


7.711 


6.217 


Ml I 


M03 


A 



Table 1 — Continued 



No. 


Name 


R.A.(J2000) 


Decl.(J2000) 


J 


H 


K S 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


Spectral Type a 


Reference 15 


Sourc 


73 


[M2002]SMC048122 


01 


00 09.42 


-72 


08 44.5 








8.800 


8.931 


8.810 


8.692 


8.200 


Kl I 


M03 


A 


74 


[M2002]SMC049033 


01 


00 30.43 


-71 


58 24.7 


9.741 


9.016 


8.775 


8.589 


8.689 


8.563 


8.476 


8.322 


K5 I 


M03 


A 


75 


[M2002]SMC049428 


01 


00 40.32 


-72 


35 58.8 


10.240 


9.511 


9.284 


9.113 


9.225 


9.077 


9.021 


8.927 


KO-7 I 


M03 


A 


76 


N371 R20 


01 


00 41.56 


-72 


10 37.0 


8.955 


8.220 


7.963 


7.821 


7.827 


7.620 


7.446 


7.207 


K5-M0 I 


WBF83 


A 


77 


HV11417 


01 


00 48.17 


-72 


51 02.1 


9.754 


8.806 


8.453 


7.891 


7.795 


7.568 


7.108 


4.686 




V08 


A 


78 


[M2002]SMC049990 


01 


00 54.13 


-72 


51 36.3 


9.517 


8.924 


8.677 


8.687 


8.682 


8.417 


8.248 


7.389 


K5 I 


M03 


A 


79 


[M2002]SMC050028 


01 


00 55.20 


-71 


37 52.9 


9.078 


8.289 


8.043 


7.458 


7.580 


7.176 


7.064 


5.406 




M07 


A 


80 


[M2002]SMC050237 


01 


01 00.31 


-72 


13 41.6 


10.413 


9.623 


9.459 


9.336 


9.440 


9.326 


9.250 


9.267 


K2-5 I 


M03 


A 


81 


[M2002]SMC050348 


01 


01 03.26 


-72 


04 39.4 


10.290 


9.583 


9.363 


9.209 


9.298 


9.188 


9.099 


8.963 


K7 I 


M03 


A 


82 


[M2002]SMC050360 


01 


01 03.58 


-72 


02 58.5 


10.427 


9.570 


9.388 


9.259 


9.391 


9.262 


9.208 


9.135 


K5 I 


M03 


A 


83 


[M2002]SMC050840 


01 


01 15.99 


-72 


13 10.0 


9.386 


8.576 


8.368 


8.118 


8.187 


8.011 


7.803 


6.131 


Ml I 


M03 


A 


84 


[M2002]SMC051000 


01 


01 19.92 


-72 


05 13.1 


10.224 


9.416 


9.237 


9.104 


9.198 


9.065 


8.981 


8.887 


K5 I 


M03 


A 


85 


[M2002]SMC051265 


01 


01 26.89 


-72 


01 41.3 


10.226 


9.414 


9.235 


9.582 


9.201 


9.057 


8.981 


8.787 


K3-5 I 


M03 


A 


86 


HV11427 


01 


01 27.80 


-72 


48 29.0 


12.660 


11.821 


11.504 


10.914 


10.747 


10.412 


10.344 






CF81 


M 


87 


HV1865 


01 


01 36.82 


-72 


44 41.1 


10.438 


9.550 


9.205 


8.877 


8.838 


8.524 


8.336 


8.147 




CF81 


M 


88 


[M2002]SMC051906 


01 


01 43.57 


-72 


38 25.1 


10.363 


9.596 


9.358 


9.203 


9.367 


9.229 


9.160 


9.076 


K5 I 


M03 


A 


89 


[M2002]SMC052334 


01 


01 54.16 


-71 


52 18.8 


9.689 


8.962 


8.695 


8.558 


8.669 


8.502 


8.365 


7.321 


K5-M0 I 


M03 


A 


90 


[M2002]SMC052389 


01 


01 55.43 


-72 


00 29.5 


9.912 


9.104 


8.910 


8.727 


8.872 


8.717 


8.604 


8.143 


K7 I 


M03 


A 


91 


[M2002]SMC053557 


01 


02 23.71 


-72 


55 21.2 


10.008 


9.241 


9.024 


8.869 


9.063 


8.882 


8.815 


8.736 


K7 I 


M03 


A 


92 


[M2002]SMC053638 


01 


02 25.83 


-72 


38 56.9 


10.398 


9.563 


9.370 


9.180 


9.311 


9.157 


9.085 


9.153 


K2-7 I 


M03 


A 


93 


[M2002]SMC054111 


01 


02 37.22 


-72 


16 25.1 


9.872 


9.149 


8.969 


8.811 


8.889 


8.773 


8.697 


8.573 


KO-5 I 


M03 


A 


94 


[M2002]SMC054300 


01 


02 42.12 


-72 


37 29.1 


10.243 


9.437 


9.224 


9.031 


9.196 


9.073 


9.007 


8.951 


KO-5 I 


M03 


A 


95 


[M2002]SMC054414 


01 


02 44.82 


-72 


01 51.9 


10.340 


9.574 


9.386 


9.217 


9.360 


9.238 


9.184 


8.740 


K3-5 I 


M03 


A 


96 


[M2002]SMC054708 


01 


02 51.37 


-72 


24 15.5 


10.055 


9.296 


9.081 


8.957 


9.055 


8.882 


8.804 


8.574 


Kl I 


M03 


A 
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No. Name R.A.(J2000) Decl.(J2000) J H K s [3.6] [4.5] [5.8] [8.0] [24] Spectral Typc a Reference b Source c 



97 


N371#29 


01 


03 


02 


.38 


-72 


01 


52. 


9 


9.901 


9.011 


8 


.619 


8. 


.328 


8 


.136 


7 


.877 


7. 


.530 


5 


.890 


M4.5 I 


WBF83 


A 


98 


[M2002JSMC055275 


01 


03 


04 


.34 


-72 


34 


12. 


.8 


9.737 


8.890 


8 


.654 


8. 


.433 


8 


.536 


8 


.311 


8. 


.167 


7 


.974 


K7-M0 I 


M03 


A 


99 


[M2002JSMC055355 


01 


03 


06 


.43 


-72 


28 


35. 


1 


9.836 


9.049 


8 


.842 


8. 


.706 


8 


.654 


8 


.490 


8. 


.289 


8 


.158 


K7-M0 I 


M03 


A 


100 


[M2002] SMC055470 


01 


03 


08 


.88 


-71 


55 


50 


.8 


10.427 


9.716 


9 


.548 


9 


.367 


9 


.436 


9 


.306 


9. 


.208 


9 


.161 


K3-5 I 


M03 


A 


101 


[M2002]SMC055560 


01 


03 


10 


.93 


-72 


18 


32 


9 


10.227 


9.444 


9 


.225 


9 


.114 


9 


.237 


9 


.069 


8. 


.987 


8 


.852 


K3-5 I 


M03 


A 


102 


N371 C12 


01 


03 


12 


.98 


-72 


09 


26 


5 


9.642 


8.852 


8 


.594 


8. 


.423 


8 


.423 


8 


.214 


7. 


.960 


6 


.132 


K5-M0 I 


WBF83 


A 


103 


[M2002] SMC055933 


01 


03 


18 


.56 


-72 


06 


46. 


2 


10.095 


9.357 


9 


.170 


8. 


.932 


8 


.982 


8 


.804 


8. 


.523 


7 


.622 


K3-5 I 


M03 


A 


104 


[M2002] SMC056389 


01 


03 


27 


.61 


-72 


52 


09 


.4 


8.802 


8.059 


7 


.779 


7. 


.580 


7 


.647 


7 


.445 


7. 


.240 


5 


.576 


K5-7 I 


M03 


A 


105 


[M2002] SMC056732 


01 


03 


34 


.30 


-72 


06 


05. 


.8 


9.940 


9.122 


8 


.864 


8. 


.683 


8 


.774 


8 


.630 


8. 


.523 


8 


.437 


K5-M0 I 


M03 


A 


106 


HV11452 


01 


03 


36 


.92 


-73 


33 


37 


.8 


10.887 


10.023 


9 


.689 


9 


.496 


9 


.480 


9 


.296 


9. 


.056 


8 


.877 




CF81 


A 


107 


[M2002] SMC057386 


01 


03 


47 


.35 


-72 


01 


16. 





10.002 


9.216 


8 


.968 


8. 


.854 


8 


.982 


8 


.862 


8. 


.776 


8 


.734 


Kl I 


M03 


A 


108 


[M2002] SMC057472 


01 


03 


48 


.89 


-72 


02 


12 


.7 


10.003 


9.207 


8 


.967 


8. 


.836 


8 


.929 


8 


.814 


8. 


.712 






K2 I 


M03 


A 


109 


[M2002JSMC058149 


01 


04 


02 


.77 


-72 


05 


27 


.7 


10.247 


9.474 


9 


.257 


9 


111 


9 


.243 


9 


.094 


9. 


.021 


8 


.813 


K2-5 I 


M03 


A 


110 


[M2002]SMC058472 


01 


04 


09 


.52 


-72 


50 


15. 


3 


10.320 


9.578 


9 


.352 


9 


.261 


9 


.249 


9 


.038 


8. 


.910 


8 


.156 


K0 I 


M03 


A 


111 


HV1956 


01 


04 


15 


.46 


-72 


45 


19 


9 


10.474 


10.069 


9 


.910 


9 


.313 


9 


.321 


9 


.174 


9 


.028 


9 


.018 


K2 I 


CF81 


A 


112 


[M2002]SMC058839 


01 


04 


17 


.71 


-71 


57 


32. 


5 


10.401 


9.625 


9 


.413 


9 


.224 


9 


.373 


9 


.229 


9 


.152 


9 


.091 


K2 I 


M03 


A 


113 


HV1963 


01 


04 


26 


.63 


-72 


34 


40. 


3 


10.638 


9.677 


9 


.409 


9 


.077 


9 


.210 


9 


.011 


8. 


.916 


8 


.571 




CF81 


A 


114 


[M2002]SMC059426 


01 


04 


30 


.26 


-72 


04 


36. 


1 


10.006 


9.174 


8 


.963 


8. 


.698 


8 


.851 


8 


.688 


8. 


.549 


8 


.053 


K5-7 I 


M03 


A 


115 


[M2002]SMC059803 


01 


04 


38 


.16 


-72 


01 


27. 


2 


9.093 


8.301 


8 


.097 


7 


.858 


7 


.860 


7 


.620 


7. 


383 


6 


.432 


K2-3 I 


M03 


A 


116 


[M2002]SMC060447 


01 


04 


53 


.05 


-72 


47 


48. 


.5 


10.266 


9.403 


9 


.175 


9 


.018 


9 


.112 


8 


.916 


8. 


.815 


7 


.860 


K2 I 


M03 


A 


117 


HV12179 


01 


05 


01 


.76 


-72 


06 


56 


.4 


10.951 


10.110 


9 


.851 


9 


.522 


9 


.409 


9 


.230 


9 


.008 


8 


.480 




WBF83 


A 


118 


[M2002]SMC061296 


01 


05 


11 


.50 


-72 


02 


27. 


5 


10.193 


9.437 


9 


.237 


9 


.061 


9 


.129 


8 


.959 


8. 


.867 


8 


.882 


K7 I 


M03 


A 


119 


[M2002]SMC062427 


01 


05 


40 


.04 


-71 


58 


46 


.4 


10.396 


9.641 


9 


.435 


9 


.288 


9 


.422 


9 


.306 


9 


.233 


9 


.203 


K5 I 


M03 


A 


120 


[M2002]SMC063114 


01 


06 


01 


.37 


-72 


52 


43 


.2 


9.800 


9.037 


8 


.788 


8. 


.579 


8 


.619 


8 


.485 


8. 


321 


7. 


.997 


K5-7 I 


M03 


A 



Table 1 — Continued 



No. 


Name 


R.A.(J2000) 


Decl.(J2000) 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


Spectral Type a 


Reference 13 


Source 


121 


[M2002]SMC063131 


01 


06 


01.72 


-72 


24 


03.8 


10.511 


9.729 


9.563 


9.367 


9.454 


9.326 


9.241 


9.110 


K5 I 


M03 


A 


122 


[M2002]SMC063188 


01 


06 


03.21 


-72 


52 


16.0 


10.202 


9.428 


9.187 


8.951 


9.050 


8.899 


8.793 


8.428 


K2 I 


M03 


A 


123 


[M2002]SMC064448 


01 


06 


40.21 


-72 


28 


45.2 


10.407 


9.770 


9.583 


9.486 


9.522 


9.416 


9.329 


9.278 


K2-7 I 


M03 


A 


124 


[M2002]SMC064663 


01 


06 


47.62 


-72 


16 


11.9 


9.241 


8.520 


8.312 


7.995 


8.033 


7.786 


7.456 


5.946 


K7 I 


M03 


A 


125 


[M2002]SMC066066 


01 


07 


29.36 


-72 


30 


45.7 


10.036 


9.266 


9.090 


8.948 


9.078 


8.920 


8.873 


8.846 


K5 I 


M03 


A 


126 


[M2002]SMC066694 


01 


07 


48.88 


-72 


23 


42.4 


9.790 


8.977 


8.774 


8.691 


8.764 


8.609 


8.581 


8.508 


K5 I 


M03 


A 


127 


[M2002]SMC067509 


01 


08 


13.34 


-72 


00 


02.9 


10.126 


9.344 


9.172 


9.047 


9.148 


9.038 


8.974 


8.918 


Kl I 


M03 


A 


128 


[M2002]SMC067554 


01 


08 


14.65 


-72 


46 


40.8 


10.129 


9.396 


9.174 


8.989 


9.126 


8.959 


8.867 


8.498 


K5-7 I 


M03 


A 


129 


[M2002]SMC068648 


01 


08 


52.08 


-72 


23 


07.0 


9.683 


8.918 


8.716 


8.579 


8.730 


8.583 


8.506 


8.420 


K2 I 


M03 


A 


130 


HV12956 


01 


09 


02.20 


-71 


24 


10.0 


11.532 


10.853 


10.336 


9.414 


9.134 


8.568 


7.301 


3.143 




CF81 


A 


131 


[M2002]SMC069886 


m 


no 
U9 


oo.Uo 


- ( O 


on 


ui.y 


8.849 


8.136 


7.793 


7.531 


7.373 


7.050 


6.758 


4.800 


K5-M0 I 


L06 


A 


L-iZ. 


UV91 1 9 


01 


10 


04.20 


-72 


36 


54.0 




















or oi 


A 


133 


HV 859 


01 


10 


27.10 


-72 


35 


48.0 


10.278 


9.522 


9.083 


8.579 


8.483 


8.258 


8.051 


7.990 




CF81 


A 


134 


[M2002]SMC071507 


01 


10 


50.25 


-72 


00 


14.5 


10.554 


9.822 


9.595 


9.481 


9.607 


9.502 


9.444 


9.242 


K3-5 I 


M03 


A 


135 


[M2002]SMC071566 


01 


10 


53.51 


-72 


25 


40.0 


10.365 


9.594 


9.393 


9.194 


9.327 


9.201 


9.152 


8.979 


K7 I 


M03 


A 


136 


[M2002]SMC081668 


01 


21 


54.03 


-73 


26 


49.2 


10.304 


9.549 


9.355 


9.133 


9.295 


9.153 


9.056 


8.883 




M03 


A 


137 


[M2002]SMC081961 


01 


25 


38.80 


-73 


21 


55.6 


9.031 


8.217 


7.984 


7.942 


8.023 


7.859 


7.744 


7.059 




M03 


A 


138 


[M2002]SMC082159 


01 


26 


09.91 


-73 


23 


15.4 


10.069 


9.351 


9.167 


9.023 


9.129 


9.007 


8.929 


8.911 




M03 


A 


139 


[M2002]SMC083202 


01 


29 


18.52 


-73 


01 


59.3 


8.861 


8.178 


7.949 


7.816 


7.894 


7.722 


7.603 


7.465 




M03 


A 


140 


[M2002]SMC083593 


01 


30 


33.92 


-73 


18 


41.9 


9.636 


8.955 


8.601 


8.020 


7.811 


7.509 


7.286 


5.706 


M2 I 


M03 


A 



a Spectral types come from Masscy & Olscn (2003); Lcvcsquc ct al. (2006, 2007) and all have been updated. 

b Reference: 'CF81' for Icatchpole fc Feastl l|l981l ). 'WBF83' for Iwood. Bessell fc Foxl jl983T ). 'M03' for llVlassev fc Olsenl J2OO3J), 'L06' for iLevesaue et al. 



(2006), 'V08' for 



van Loon et al 



(2008), 'V10' for 



van Loon et al 



(2010). 



c Photomctry data resource: 'A' for ASAS, 'M' for MACHO and ' — ' for null. 
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Table 2. Results of robust linear fitting for relationship between infrared and bolometric 

magnitude 



Band 


Slope a 


Intercept b 


a 


J 


0.65±0.05 


2.82±0.17 


0.31 


H 


0.63±0.05 


2.23±0.18 


0.33 


K S 


0.65±0.06 


1.78±0.19 


0.35 


[3.6] 


0.69±0.06 


1.23±0.21 


0.38 


[4.5] 


0.72±0.06 


0.97±0.21 


0.37 


[5.8] 


0.79±0.07 


0.07±0.22 


0.40 


[8.0] 


0.86±0.07 


-0.86±0.24 


0.44 


[24] 


1.45±0.14 


-7.82±0.48 


0.86 



Table 3. Results of time-series analysis by PDM, Period04 and CLEAN methods for 6 semi-regular and 15 LSP 

RSGs 



Method PDM Period04 CLEAN 



No. 


(i»v) 


M b oi 


P (d) 


A (mag) 





P (d) 


A (mag) 


Phase 


P (d) 


<P>±<x(d) 


Points 


Type a Ref b 


5 


13.11 


-8.48 


600 


0.33 


0.35 


620 


0.22 


0.84 


623 


614±13 


208 


L 


L06 


8 


13.01 


-8.38 


593 


0.31 


0.62 


603 


0.16 


0.13 


595 


597±5 


208 


L 


L06 


10 


12.11 


-8.47 


572 


0.17 


0.41 


566 


0.15 


0.01 


569 


569±3 


214 


S 


L06 


12 


12.51 


-7.93 


366 


0.17 


0.46 


358 


0.08 


0.14 


363 


362±4 


214 


L 


L06 


21 


12.86 


-8.13 


388 


0.18 


0.55 


381 


0.08 


0.75 


385 


384±4 


212 


S 


L06 


29 


12.82 


— 


512 


0.30 


0.43 


510 


0.31 


0.66 


520 


514±5 


196 


S 




30 


12.42 


-8.39 


568 


0.23 


0.60 


562 


0.09 


0.97 


569 


566±4 


205 


L 


L06 


37 


12.61 


— 


620 


0.17 


0.38 


618 


0.20 


0.32 


623 


620±3 


205 


S 




41 


11.96 


-8.44 


531 


0.14 


0.44 


513 


0.13 


0.78 


523 


522±9 


197 


L 


L06 


42 


12.12 


— 


483 


0.13 


0.30 


484 


0.07 


0.19 


484 


483±1 


206 


L 




50 


12.57 




309 


0.11 


0.61 


313 


0.07 


0.08 


319 


313±5 


200 


S 




54 


13.05 


-7.88 


419 


0.24 


0.40 


427 


0.16 


0.73 


434 


426±8 


186 


L 


L06 


56 


12.07 




195 


0.11 


0.67 


195 


0.04 


0.17 


195 


195±0 


195 


S 




72 


12.51 




517 


0.29 


0.46 


522 


0.23 


0.11 


542 


527±13 


189 


L 




79 


12.36 


-8.80 


721 


0.44 


0.30 


716 


0.42 


0.07 


723 


720±4 


190 


L 


L09 


83 


12.61 


-8.12 


541 


0.16 


0.45 


571 


0.11 


0.48 


566 


559±16 


187 


L 


L06 


98 


12.47 




464 


0.21 


0.49 


448 


0.08 


0.19 


448 


453±9 


182 


L 




111 


12.35 




210 


0.34 


0.17 


209 


0.29 


0.66 


210 


209±1 


189 


L 




115 


12.17 


-8.69 


541 


0.22 


0.33 


534 


0.13 


0.22 


566 


547±17 


189 


L 


L06 


131 


12.08 


-8.76 


728 


0.39 


0.26 


728 


0.36 


0.08 


766 


740±22 


190 


L 


L06 


140 


12.54 


-8.06 


506 


0.50 


0.46 


499 


0.29 


0.12 


517 


507±9 


175 


L 


L07 



a Type: 'S' for semi-regular RSGs, 'L' for LSP RSGs with distinguish short period 
b 'Ref of M bol : 'L06' for 



Levesque et al 



(2006), 'L07' for 



Levesque et al 



(2007), 'L09' for 



Levesque et al 



(2009). 
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Table 4. K s -band P-L relations of RSGs in the Galaxy, LMC, SMC and M33 



Galaxies 


Slope a 


intercept b 


a 


Galaxy 


-4.34±0.60 


0.91±0.98 


0.51 


LMC 


-3.74±0.30 


-0.38±0.51 


0.24 


SMC 


-3.28±0.39 


-1.66±0.64 


0.27 


M33 


-3.59±0.41 


-0.88±0.69 


0.31 


Combination 


-3.68±0.18 


-0.62±0.30 


0.31 
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Table 5. Near- and Mid-infrared P-L relations of RSGs in SMC and LMC 



Band 




SMC 






LMC 




SMC+LMC 




Slope a 


intercept b 


a 


Slope a 


intercept b 


a 


Slope a 


intercept b 


a 


J 


-2.57±0.35 


-2.56±0.58 


0.25 


-3.61±0.38 


0.39±0.63 


0.30 


-3.15±0.25 


-0.93±0.42 


0.28 


H 


-3.08±0.39 


-1.96±0.63 


0.27 


-3.58±0.33 


-0.50±0.54 


0.26 


-3.30±0.23 


-1.30±0.38 


0.26 


K S 


-3.28±0.39 


-1.66±0.64 


0.27 


-3.75±0.30 


-0.38±0.51 


0.24 


-3.55±0.22 


-0.90±0.37 


0.25 


[3.6] 


-3.11±0.37 


-2.37±0.60 


0.25 


-3.97±0.28 


-0.12±0.47 


0.22 


-3.78±0.22 


-0.64±0.36 


0.24 


[4.5] 


-3.29±0.35 


-1.87±0.57 


0.24 


-4.34±0.25 


0.89±0.41 


0.20 


-4.10±0.20 


0.26±0.33 


0.22 


[5.8] 


-3.58±0.37 


-1.33±0.61 


0.26 


-4.53±0.28 


1.13±0.46 


0.22 


-4.41±0.22 


0.83±0.36 


0.25 


[8.0] 


-4.03±0.39 


-0.33±0.64 


0.27 


-5.33±0.39 


2.90±0.64 


0.31 


-5.36±0.30 


3.07±0.49 


0.30 


[24] 


-7.73±1.04 


8.50±1.70 


0.72 


-7.82±0.69 


7.80±1.16 


0.56 


-9.36±0.66 


12.34±1.10 


0.74 
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Table 6. Period of LSP of 31 SMC RSGs obtained by PDM, Period04 and CLEAN 



No. A (mag) PDM (d) Period04 (d) CLEAN (d) (LSP) ± a (d) 



5 


0.19 


2460 


2436 


2380 


2425+41 


6 


0.06 


1828 


1816 


1870 


1838+28 


8 


0.19 


2406 


2445 


2510 


2453+53 


12 


0.07 


1588 


1553 


1625 


1 KOQ-LOC 

-L ' lUO 1 -Ml 


14 


0.14 


2697 


2656 


2598 


2650+50 


24 


0.02 


1720 


1549 


1613 




26 


0.17 


1627 


1606 


1636 


i R23+1 5 


32 


0.04 


2024 


2017 


2182 


2074+93 


42 


0.13 


1848 


1811 


1793 


1817+28 


45 


0.11 


2582 


2474 


2633 


' H >■ > 1 (1 J_ 


54 


0.11 


1283 


1295 


1302 


i 2QS+1 


62 


0.08 


2589 


2545 


2566 


2566+22 


64 


0.07 


1615 


1601 


1624 


1613+12 


69 


0.10 


2549 


2538 


2507 


2531+22 


74 


0.04 


2010 


2045 


2002 


2019+23 


79 


0.36 


2089 


2007 


2145 


2080+69 


83 


0.11 


1849 


2004 


2151 


2001+151 


84 


0.07 


2749 


2521 


2604 


2624+115 


89 


0.09 


1620 


1610 


1640 


1623+15 


93 


0.04 


1420 


1309 


1434 


1387+68 


95 


0.09 


2533 


2546 


2582 


2553+25 


98 


0.08 


2212 


2294 


2230 


2245+43 


100 


0.06 


1385 


1377 


1447 


1403+38 


107 


0.04 


1702 


1740 


1807 


1749+53 


110 


0.13 


1482 


1574 


1613 


1556+67 


114 


0.04 


1827 


1891 


1845 


1854+33 


115 


0.10 


2507 


2540 


2582 


2543+38 


119 


0.07 


2414 


2537 


2575 


2508+84 


120 


0.05 


1154 


1209 


1184 


1182+28 


128 


0.04 


1849 


1746 


1860 


1818+63 


129 


0.03 


1620 


1600 


1628 


1616+14 



- 51 - 



-52 - 



-69.000 



-70.000 



-71.000 



-72.000 



-73.000 



-74.000 



-75.000 



-76.000 



-77.000 




40.000 30.000 20.000 

Right Ascension (degree) 



10.000 
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Fig. 1. — The spatial distribution of all the 140 sample stars (red circles) superposed on the 
Spitzer /SAGE 8 /zm mosaic image, mostly in the bar region. 
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Fig. 2.— The relationsh ip between the infrare d and bolometric mag nitud es of the 113 



reliab le RSGs targets from 



Massev fc Olsenl (120031 ) 



Levesque et al. 



(j2006f ) and 



Massey et al. 



(120071 ). The red solid line shows the robust linear fitting. The re sult in the Ks band indicates 



a moderate discrepancy between ours and 



Josselin et al. 



((2000|)'s (black dash-line) possibly 



due to the influence of metallicity and extinction. 
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Fig. 3. — The Kg vs . J — K g CMD for all the targets. For comparison, the 3654 massive stars 



from 



Bonanos et al 



( I2010I ) are added as background and denoted by black dots. Different 
symbols denote different resources. Most of the targets have occupied the region of 7.5 
< Ks < 10 mag and 0.5 < J — Ks < 1.6 mag. The vertic al dash-line at J — Kq = 1.6 



mag indicates the boundary of carbon-rich stars suggested by 



Wood. Bessell & Foxl (Jl983|) . 



The other dash-line at J — Kg = 0.5 mag shows the observational boundary of RSGs by 



Josselin et al. 



( l2000l ). The horizontal dash-lines show the theoretical luminosity boundaries 



in Kg band with defined masses. The vertical dot-line indic a tes th e limit of evolutionary 



Origlia et al. 



( 119991 ). The horizontal dot- 



synthesis models of RSGs in near-infrared by 
line indicates the cri t erion of Mboi = —7.1 mag used to distinguish AGBs and RSGs by 



Wood. Bessell fc Foxl ( 119831 ). Obvious outliers are marked by their ID in Table [TJ The 



conventions of the symbols for resources are kept in following figures. 
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Fig. 4. — The same as Fig. |3] 



Dut for different classifications of evolved stars. The J — Kq cuts 



shown as solid lines are from 



Cioni et al. 



fl2006a 



bj), where C-AGBs are defined by K s < KO 

re red giant 



and K s > K2, O-AGBs by K s < KO and Kl < K s < K2. The K s -band tip of t 
branch (TRGB) is denoted by dot-line. The profiles of contour maps are from 



Boyer et al. 



( 120 111 ) with red, blue and purple solid lines for RSG, O-AGBs and C-AGBs respectively. 



For a better view of distinguishing RSGs from AGBs, all the O-A GBs (dark b l ue cir cle) 



Boyer et al. 



fboilh are 



C-AGBs (purple circle) and x-AGBs (green circle) from Table 4 of 
overplotted. From the diagram, the RSGs sample of Boyer has too low luminosity and too 
narrow and blue color index, where the red dot-line shows the region of RSGs after moving 
redward 0.1 mag and upward 1.5 mag. 
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Fig. 5. — T he [3.6] vs. 



and RSGs by 



Boyer et al. 



Bonanos et al. 



(j20inh 



— [3.6 ] diagram. The regions of AGBs by 

(1201 ll ) are denoted by black and red solid lines respectively. The 
limits of [3.6]-band magnitude are shown as horizontal dash-lines according to the RSGs 
masses. The boundaries of color index are J — [3.6 ] = 0. 59 mag and J — [3.6] = 2.62 mag, 



i.e. the bluest point o f RSGs region by 



Bonanos et all (I2010h . 



Boyer et al. 



( 1201 ll ) and the reddest edge of AGBs by 
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Fig. 6.— The [8.0] vs. K s — [8.0] diagram. The [8.0] magnitude limits are set up at 
[8.0] = 6.78 mag (30 M ) and [8.0] = 9.90 mag (13 M ). The K s - [8.0] gets larger with 
increased [8.0], which can be fitted by a line as [8.0] = (-1.05±0.11)(K s -[8.0]) + (9.08±0.10) 
with a = 0.71 shown as a dot-line. According to this line, the red limit of Kg — [8.0] is decided 
as the cross-point of the line with the upper limit of luminosity. The blue limit of Kg — [8.0] 
is set at 0.1 mag since no targets bluer than it. 
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Fig. 7.— The [8.0] vs. J - [8.0] diagram. As Fig. El the linear fit [8.0] = (-1.0 ± 0.88)(J - 
[8.0]) + (10.08 ± 0.13) with a = 0.70 is shown as a dot-line. The limits of J - [8.0] are set at 
3.31 and 1.0 mag by the same way as in Fig. |6j The dashed line at [8.0 ] = 8.4 mag is the 
upper limit of AGBs in the [8.0]-band (see Figure 4 in 



Boyer et al. 



2011( 1 . 
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Fig. 8. — The same as Fig. [61 but for [24] vs. [8.0] - [24]. The correlation of the brightness 
in [24] with the color index [8.0] — [24] is very tight since both reflect the radiation of the 
dust. The linear fit for this relationship is [24] = (-1.96 ± 0.05)([8.0] - [24]) + (9.03 ± 0.06) 
with a = 0.45. 
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Fig. 9. — The two-color diagram J — Ks versus Ks — [8.0]. The black solid and dash-line 
represents the black bod y (BB) model at diffe rent temperatures and a mod el of a 3500 K 



BB p lus 250 K dust from 



Bonanos et al 



Bonanos et al. 



2009 



(|2010| ) respectively (see Section 4 of 
for detailed description). The limits of 0.5 ~ 1.6 mag in J — Ks and 0.1 ~ 1.17 mag 
in K s — [8.0] are shown by dot-lines respectively. The outliers are the same as in Fig. [6] and 
not marked. 
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Fig. 10. — The same as Fig. HI but for [3.6] — [4.5] vs. [4.5] — [8.0]. Due to the impact of 
the CO line in the [4.5] band, most targets are far from the BB model. There is no certain 
limits for these color indexes, and only the obvious outliers are labeled by their IDs. 
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Fig. 11. — The same as Fig. [9j but for [3.6] — [8.0] vs. [8.0] — [24]. As the circumstellar 
dust dims the [8.0]-band brightness, most targets are below the BB line but approximately 
follow the 3500 K BB plus 250 K dust model. The outliers in previous diagrams also roughly 
follow this model. Only Nos.4 and 130 are far from the model and marked. 
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Fig. 12.— Histogram of J - K s of RSGs in SMC (solid line) and LMC (dash-line). The 
difference of J — K s between them is about 0.1 mag. 
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Fig. 13. — The CMD for absolute magnitude in the Ks-band versus J — Ks of RSGs in SMC 
(red triangle) and LMC (black diamond). The J — Ks in SMC is shifted redward by 0.1 mag 
to compensate for the metallicity effect. The SMC RSGs are slightly fainter and bluer than 
the LMC. 
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Fig. 14. — An example of processing the ASAS light curve. Top to bottom: original light 
curve; least-square (Savitzky-Golay) polynomial smoothed light curve; 10 day-binned light 
curve. 
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Fig. 15. — The MACHO light curves. Three outliers that are too faint to be RSGs in the 
left column and three other similar targets in the right column are illustrated and shown 
with their ID number, period and amplitude of variation. 
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Fig. 16. — The process to derive the period using the PDM2 method based on the ASAS data 
for the object No. 12 which used also in following diagrams. Top to bottom: Left column : 
10 day-binned light curve, theta diagram, phase diagram; right column: residual light curve, 
residual theta diagram, residual phase diagram. The dash-lines in the theta diagrams are 
significance levels of 0.05 and 0.01 respectively. The red lines in the bottom panel are the 
fitted curve. 
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Fig. 17. — The same as Fig. [161 but using the Period04 method. Top panel is the 10 day- 
binned light curve with the fitted curve; bottom panels are power spectra. The spectral 
window in the first power spectrum is shown at the same scale as the power spectra. The 
dash-line shows the 4<r level. The red cross marks the highest peak in the power spectra in 
each iteration. 
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Fig. 18.— An example of the CLEAN processing. From top to bottom: the "dirty" 
spectrum, clean components, residual spectrum, clean spectrum. The red cross marks the 
highest peak in the clean spectrum. 
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Fig. 19. — An example of the WWZ processing. From top to bottom: the contour map 
of WWZ, the time-varying amplitude, the time-varying period. The dense regions in the 
contour map are marked by dot, dash and solid lines at the level of 32, 64 and 128 respectively. 
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Fig. 20. — The typical light curves of three kinds of variation in the ASAS data and two 
RSGs with the MACHO data. Top panel from left to right: irregular, semi-regular and 
long secondary period variables. Bottom panel: two MACHO targets with complex light 
variation. In each diagram, the target ID is present in the top left corner, and red line for 
the fitted curve. 



-72 - 



0.5 F 




Log P 

Fig. 21. — The change of amplitude with period for the 6 semi-regular RSGs and 15 LSP 
RSGs with distinguishable short period, where the red solid line is the linearly fitted curve, 
AV = (0.76 ± 0.14) x logP - (1.83 ± 0.24) with a = 0.10. 
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Fig. 22. — The Kg-band P-L relation of the RSGs in SMC, where the fitted line (red solid) 
writes K s = (-3.28 ± 0.39) x logP + (17.25 ± 0.64) with a = 0.27. 
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Fig. 23. — The same as Fig. [22] but for absolute magnitude in the Ks-band and the RSGs 
in the Galaxy, LMC, SMC and M33 respectively. The convention of symbols and line styles 
are kept in following diagrams. The fitted parameters are listed in Table. H] as a form of 
K s = a x log P + b. 
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Fig. 24.- The K s -band P-L relation of the RSGs in the Galaxy, LMC, SMC and 



M33, similar 
other people ( 


;o Fie. I231 


with 


the addition of the P- 


j relations of RSGs derived by 


Feast et al. 




L980 




CatchDole & Feast 


198] 




Kinman. Mould & Wood 


1987; 


Pierce. Jurcevic & Crabtree 


2000; 


Kiss. Szabo & Beddine 


2006 


) . The AGB 02 sequence from 


Soszvhski et al 


( 


2007 


) is almost superposed on ours. 
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Fig. 25. — Comparison of the P-L relation in the V, J, H, and Kg bands of RSGs with the 
LMC The black dash- and dot-lines are the fitted line for SMC and LMC respectively, and 
the red solid line for a combination of SMC and LMC. The V band radiation suffers the 
interstellar and circumstellar extinction and exhibits no correlation with the period. 
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Fig. 27. — The P-L relations for the LSP of LSP RSGs where the short periods are plotted 



for comparison. 



Soszyhski et al 



Super posed in this diagram are the P-L relati ons of LPVs (so 



( 120071 ) and various LPVs targets (points) from 



Soszyhski et al 



i d line s) from 
(ho 111 ) in the 



SMC, where the sym bols are the sam e as t 



and fainter objects in 



Soszyhski et al. 



Soszyhski et al. 



re original papers. With the addition of more 



(I2011I ). there are small shifts of the P-L relations by 



(120071 ) from the locations of the objects. 
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Fig. 28. — The P-L relation for the first three radial modes (fundamental mode at the right 
an d second o v erton e at the left) at 15, 20, and 30 M stars with Z=0.005 from the model 



by 



Guo fc Lil (120021 ) . The objects in our sample with reliable luminosity are shown mostly 



between 20 and 30 M Q by squares. 
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Fig. 29. — The same as Fig. [29] but for the models with Z=0.01 from 
RSGs in LMC. 
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Fig. 30.— The same as Fig. [29] but for the models at Z=0.02 from lOuo fc Lil (J2002) and 
RSGs in the Galaxy. 
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Fig. 31. — The Same as Fig. 1291 but for the models by iHeger et al.l (119971 1 which have no 
discrimination of the metallicity and for all the objects in SMC, LMC and the Galaxy. 



